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ABSTRACT 


The rate of convergence of infinite series can be accelerated by a suitable 
splitting of each term into two parts and then combining the second part 
of the n-th term with the first part of the (n+l)-th term to get a new series 
and leaving the first part of the first term as an "orphan”. Repeating this 
process an infinite number of times, the series will often approach zero, and 
we obtain the series of orphans, which may converge faster than the original 
series. Heuristics for determining the splits are given. Various mathematical 
constants, originally defined as series having a term ratio which approaches 
1, are accelerated into series having a term ratio less than 1. This is done 
with the constants of Euler and Catalan. The series for pi/4 = arctan 1 is 
transformed into a variety of series, among which is one having a term ratio 
of 1/27 and another having a term ratio of 54/3125. A series for 1/pi is 
found.which has a term ratio of 1/64 and each term of which is an integer 
divided by a power of 2, thus making it easy to evaluate the sum in binary 
arithmetic. We express zeta(3) in terms of pi+3 and a series having a term 
ratio of 1/16. Various hypergeometric function identities are found, as well 
as a series for (arcsin y)+2 curiously related to a series for y arcsin y. 
Convergence can also be accelerated for finite sums, as is shown for the har¬ 
monic numbers. The sum of the reciprocals of the Fibonacci numbers has been 
expressed as a series having the convergence rate of a theta function. Finally, 
it is shown that a series whose n-th term ratio is (n+p)(n+q)/(n+r)(n+s), where 
p, q, r, s are integers, is equal to c + d pi+2, where c and d are rational. 

Work reported herein was conducted at the Artificial Intelligence Laboratory, 
a Massachusetts Institute of Technology research program supported in part 
by the Advanced Research Projects Agency of the Department of Defense and 
monitored by the Office of Naval Research under Contract Number N00014-70-A- 
0362-0005. 
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purpose of the United States Government. 
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Introduction 

This paper describes a way to manipulate sums to 

produce new ones which converge faster. For example, knowing. 

only that 

pi/4 - 1 - 1/3 + 1/5 - 1/7 + ... 

we can find that 

P> 2 112 134 1 5 G 

— “ — CB-{13 -(20 -(27 - ... 

4 . 15 3 7 11 3 13 17 3 19 23 

without even using calculus. Now the first series is a lot 
simpler, but it would require billions of terms to get 10 digits, 
of pi, whereas the second series only needs 7 terms, since 
each term is less than l/27th of the previous one. 
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Another example is the sum of the reciprocals of the cubes: 


1 1 1 

zeta 3 - — + — + — + ... 

3 3 3 

1 2 3 

Ue will shou that this s I ou series 


5 1 111 

4 1 3 4 2 


Ev 

n>l n 


2 2 1 3 3 

- - ( - - 
5 S 3 7 8 


5 

2 


E 


R-l 

{-) 


3 2k 
k ( ) 

k 


where each term is less than 1/4 of the previous. 


The treatment ui I I be informal with feu, if any, proofs. 

The method consists of splitting each series term into tuo 
pieces, then recombining the pieces of adjacent terms to 
form a neu series with smaller terms, and a leftover quantity 
called an "orphan". This process is repeated until a formula 
for the orphans is deduced, wherupon the shrinking series is 
discarded in favor of the series of orphans. The simplest 
case of this is an old trick known as 


Euler’s Transform 


Although it is usually applied to alternating series, Euler’s 
transform is more neatly derived for the general series 



a _^ 

■ 0 + > U a 

2 n 
n>0 
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a + a 
n n+1 . 

where U is the averaging operator, U a - _ 

n 2 

Thus we get the orphan term (a /2) and a new series upon which 

0 

we may repeat this operation to get 


a U a 

0 + 0 + 


z>\ 


n>0 


a + 2 a + a 

2 n n+1 n+2 . 

where U a - U (U a ) - - 

n n 4 


If we do this k times we will get 


a U a 

0 0 

2 2 


k-1 
U a 


E' 

n>0 


Now as k approaches infinity, the sum on the right approaches 0 
for all convergent (and many divergent) series, leaving, finally, 


n>0 k>0 


0 


(U a “ a ). The motivation for this transform shoulc 


become clear with an example. Ue have 
pi/4 - arctan 1 - 1 - 1/3 + 1/5 - 1/7 + ... 
n 

(-) 


so that a 


2n+l 
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n n+1 n 

1 (-) (-) (-) 

U a - - {-+-) „ - 

n 2 2n+l 2n+3 (2n+l)(2n+3) 


2 

U a 


n n+1 

1 (-) (-) 

- (-+- 

2 (2n+l)(2n+3) (2n+3)(2n+5) 


2 (-) 

(2n+l)(2n+3)(2n+5) 


It is easy to induce that 

r 

k! (-) 


k 

U a 


I2n+1) ... (2n+2k+l) 


Setting n to 0, ue find 


P' 

4 


n 

E f—) 1 V 112 3 k 

2n+l 2^2 5 7 2k+l 


n>0 k>0 

The left hand sum has as a ratio of consecutive terms 
' 2 n+1 


2n+3 


which approaches -1 as n grous large. 


To get S digits of pi uould require about half a billion terms. 
The right hand sum is easier to compute than it looks, since 


k 

its kth term is just - times the previous. 

2k+l 

Since this term ratio approaches 1/2 as k grows, only 38 terms 
are needed for S digits (30 bits). This convergence rate of 
1 bit/term is typical of Euler transformed series. 

Note, however, that if by ugly coincidence some 

different series happened to match our pi/4 series for the 

first.38 terms (and then, for instance, became all zeros), the 

first 30 orphans would have the same numerical 

values, albeit from a different formula, in which 

case few, if any, of the 9 digits would be right. 

Moral: virtually nothing is certain until you have the formula 
for the kth orphan. 
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Generalization of the Euler Transform 

Unfortunately, if the original series converges more rapidly 
than 1 bit/term (term ratio < 1/2), the transformed series 
may converge less rapidly, as ue shall see. This is a symptom 
of the fact that Euler’s transform is most effective when 
consecutive terms are of nearly equal magnitude and opposite 
sign, so that their averages will be very small, leaving most 
of the sum’s value in the orphan (More orphan than not). 

What we need is an unsymmetricaI, weighted averaging operator 
to achieve this near telescopy in series where each term is 
significantly smaller than its predecessor. We are indebted 
to Rich SchroeppeI for suggesting the weights 1/(1—r) and 
r/(r-l) be used when the series has a limiting term ratio r. 
This will actuaI Iy split the terms of non-aIternating series 
(positive r) into pieces with opposite sign. It was an easy 
step from this to realize that the weights need not be the 
same for each term. 

Thus we will denote a pair of weights by s and 1 - s , 

n n 

and call s a splitting function, 
n 

Ue generalize Euler’s transform as follows: 


E a - \ (s + 1 - s ) a 

n / -J n n n 

n>0 n>0 


■ s a + 'V ^ (1 - s ) a + s a 
0 0 / -* n n n+1 n+1 

n>0 

which specializes to the first step of Euler’s transform when 
s - 1/2. s - 0 leaves the series alone, while s =1 
n n n 

translates the series by one term, the zeroth becoming the 
orphan. Ue wish to choose s , the splitting function, so 

n 

as to make the 

summand both small and simple; small so that upon repeated 
applications of this operation, the value of the sums will 
accrue rapidly in the orphans, and simple so that these 
orphans have a general formula easy to discover and compute. 
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For examples, let s = s, a constant. Transforming 
n 

the geometric series: 


- E""" ■ * £ 


X m S 

n>0 n>0 


n n+1 

(1-s) x + s x 


s + (1 - s + s x) 


E- 


n>0 


+ (l-s + sx)s+(l-s + sx) s 


■E 


k>0 


(1 - s + s x) 


1 - (1 - s + s x) 


an identity, except the interval of convergence has moved 
from -1 < x < 1 to 1 - 2/s <x<l or 1 < x < 1 - 2/s, 
depending on the sign of s. 

I f we choose s so that 1 - s + s x => 0, all of the terms of this last 
sum vanish except the first, leaving s which is indeed l/(l-x). 

Ue note that if s = 1/2 (Euler’s transform) and 
-1/3 < x < 1/2, the transf ormed series converges more 
slowly than the original geometric series. 


Integrating these series from 0 to x; 


n k k 

- In (1 - x) - 

n>l k>l 

■ In (s - s x) - In s 
Here, if we again choose s - l/(l-x) 


E- 


s + s x) 
k 


(1 - s) 
k 
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n x k 



n>l k>l k 


or - In (1 - x) - In (1-) 

x - 1 

^ * -1/2, ue can compute In (3/2) with an asymptotic 

term ratio of 1/3 using the transformed series, versus -1/2 
for the original. 

If we had chosen s = 2/(2-x) instead of l/(l-x), then 
1-s+sx - -(1- s), that is, alternate terms in the 

sum for In (s - s x) - In s will cancel, leaving 


n x 2 k + 1 



n>l k>0 2 k + 1 


1 + y x 

or - In (1 - x) » In- with y -- 

1 - y 2 - x 


Here, if x . -1, the term ratio approaches 1/9, yielding In 2 
with almost 1 decimal digit/term instead of 1 bit/term, as 
with Euler’s transform. 


The formulas we have seen up to now have been known for a 
ong time. Before ue get into the more interesting goodies, 

I want to replace the summation notation with one that stresses 
term ratios, and reformulate the series transform accordingly. 


R Notation 

Let R r stand for the series whose first term is 1, second 
n>m n 

term is r and whose ratio of r.+lst to nth term is r 
m 

i-e., 1 + r (1 + r (1 + ... . 

m m+1 


a + a 
m m+1 


a 

a R _n+l . 
m n>m a 


Thus 










PAGE 8 


Note that m need not be an integer: 
P' 2 n + 1 


R 

n>0 2 n + 3 


1 n + 1 
n>- 
- 2 


Other examples of R notation: 

n 

x x 

6 -Enl 


n>0 


In (1 + x) 


R - 
n>l n 


£ • 


x 2 k + 1 


k>0 2 k + 1 

2 x 2 k + 1 x 2 

x + 2 k>0 2 k + 3 x + 2 


F(a+l, b+1:c+1;z) 


c! 

a! b! 


£ 


n>0 


(n+a)! <n+b) ! n 

- z 

(n+c)! n! 


n + a n + b 

■ R - - z 

n>l n+c n 

(F is the hypergeometric function. Unfortunately, a. b, and c are 

factoria'f f’ a P? arent 'y to'accomodate the Gamma notation for the 
factor a function. The author finds this unnatural and pedantic.) 

uithV^A 1 ? 5 ' 7° tation cannot express a series starting 

mth 0. Also, the (2 k + 1) / (2 k + 3) in the logarithm 
formula is awkward suggesting the following generalization: 

Let b R r stand for the series 
n n>0 n 


b + r (b 
0 0 1 


r (b 
1 2 


b +rb+rrb+rrrb+... 

0 0 1 0 1 2 0 1 2 3 
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b 

n+1 

b R .-r 

0 n>0 b n 
n 


Another way to say it is that 


b R r = S 

n n>m n m 

where S satisfies the first order, linear recurrence relation 

S * b + r S 

m mm m+1 


Then our logarithm formula becomes 


2 x x 2 

In (1 + x) - -(2 k + 1) R {-) . 

x + 2 k>0 x + 2 

Now our notation is no longer uniquely determined by 
the series—in fact one could cop out and write 

a+a + ...» a R 1 . 
m m+1 n n>m 

The intended convention, however, is that only factors of 
the nth, but not the n+lst term appear to the left of the R. 
(And, of course, the necessary scale factor to make the first 
term come out right.) This will prevent factors of the form 
(k + a)/(k + a + 1) or its reciprocal from needlessly 
complicating the expression to the right of the R, whose 
primary purpose is to get rid of nth powers and factorials. 

Thus, Ramanujan s amazing pi formula reads 

21460 n + 1123 (4 n + 1) (2 n + 1) (4 n + 3) 

- R - 

882 n >0 4 3 

128 21 (n + 1) 

(Collected Papers, US, .eq 39). 


4 

Pi 
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The advantages of this notation are; 

One can tel I at a glance the first term 
and term ratio, which immediately indicates 
the required computation effort and convergence 
rate. The number of digits/term is just 
-log Lim Ir I as n grow large, where the base 
n n 

of the log is the radix of the digits. 

The notation is usually more concise, since term 
ratios are usually simpler than terms. This often 
makes it easier to notice when a series is equal 
to, or a special case of, another. 

It suggests a generally more efficient method of 
numerically evaluating n terms of a sum— 
instead of computing the terms and adding them up, 
take r , add b , multiply by r 

n n n-1 

add b , ... multiply by r , and finally add b . 

n-1 0 0 

It has suggested tuo “spigot" methods of 
evaluating sums to unlimited accuracy without multiplying 
or dividing very long numbers. A "spigot method", is 
one that produces digits on demand, computing 

only as much as necessary, yet able to supply subsequent digits 
upon subsequent demand, without, in effect, starting over. 

It is immediately rewritable as a continued fraction: 

P 

b R _n * b 
n n>0 q 0 

n 


Although it isn t quite 3s manipulable as summation 
notation, ue can do things like 

a-1 b 

(k b + a) x R r x 
k>0 k 

Finally, It simplifies the forthcoming discussion of series 
transformat ions. 


da b 

-- x R r x 
dx k>0 k 


P 

_0 

0 + q 
_0 

b + p 
1 _1 

0 + q 


1 
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Derivation of the Transform 

After applying k-1 splitting functions, thereby extracting k-1 
orphans, we are ready to apply the kth splitting function s 

k, n 

to the remaining series, whose nth term we denote a , so that 

k, n 

a - a , the nth term of the original series. Also, 

0, n n 


a 

k, n+1 

r - - 

k,n a 

k, n 

the nth term ratio of the ktuply accelerated series. Then, 


EV 

n>m 


E a (s + 1 - s ) 
k,n k,n k,n 


- s a 
k,m k.,m 


* E“ 


s ) a 4- s a 

k,n k,n k,n+l k,n+l 


- s a 
k,m k,m 


E 

"n>m 


(1 - s + s 

k,n k,n+l 


r ) 
k, n 


a 

k, n 


and now if we define u » 1 
k, n 


s + r s 
k,n k,n k,n+l 



rt>m 


s a + V-' u a 
k»m k, m / j k, n k, n 
n>m 


so s a 
k, m k, m 


the kth orphan, and the nth term of our new series i 


u a , 
k,n k,n 


which we define to be a 

k+l,n 


Thus 


EV 

n>m 


s a + a 

k,m k,m / > k+l,n 
n>m 
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As before, we iterate on k indefinitely, choosing s so 

k,n 

that the u and thus the a and thus their sums tend to 
k,n k,n 

0 with k as rapidly as possible. Then nothing will remain of 
the original sum but the aggregate of the orphans. If k is 
i n i t i a I I y j, 


n>m 


* 8 a +s u a + s u u a 

j» m J* m j+l,m j,m j,m j+2,m j+l,m j,m j,m 

because a ■ u a by definition. 

k+l,n k,n k,n 

In R notation, this last equation reduces to 

a . R r - a s R u 

j,m n>m j,n j,m k,m k>j k,m 

Note that we have nowhere used the fact that j is initially 0, or 
even an integer, so that we are free to choose any convenient 
value to denote our original series. 

He now repeat this whole derivation in R notation, skipping fewer 
steps. The reader is advised to understand it, since the main 
results in this paper depend heavily on it. 

For convenience, we assume that the original series 

Then 


(1 - s 

j * n 


+ s R 
j,n n>m 

+ r s R r 

j,m j,n+1 n>m j,n+1 


s ) R r 
j,n n>m j, n 


j.n 


nas oeen scaled so that a = 1. 

j,m 


Z—/ J»n 


R r 
n>m j,n 


(1 - s } R r 

j,n n>m j,n 


(1 - s ) R r + s 

j,n n>m j,n j, m 


j,m j,m j+l,m j+l t m 
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(1 - s ) R r 

j.n n>m j,n 


j.n 

s +. r s R r -- 

j.m j * rn j,n+l n>m j, n+1 r 

j.n 


(1 


j.n 

s ) R r + .6 + r -s R r 

j»n n> r ii j * n j,m j, m r j, n+1 n>m 

j.m 


- (1 


s ) R r 
j,n n>m j,n 


B 

j.m 


s R r 

j.n j, n+1 n>m j,n 


s + (1 - s + r s ) R r 

j»m j.n j.n j,n+1 n>m j,n 

s + u R r 
j. m j, n n>m j, n 


u 

j. n+1 

-s + U R - r 

j.ro j.m n>m u j,n 

j.n 


- s + u R r 

j.m j.m n>m j+l,n 


k,n+l 

■ -r because a 

k+1 - n U k.n k+1, n 

k.n 


u a ) 
k.n k.n 


+ u (s + u R r ) 

j» m j* m j+l.m j +1,m n>m j.n+2 

+ u + u . (s + . 

j* m j.m j+l.m j+1,m j+2,m 


R u 

k.m k>j ' k.m 


» s 
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Final Iy, we have the 



Our problem then boils down to finding a pair of functions 
s and r which satsify the recurrence 
k, n k, n 

u 

k,n+l 

r - -r 

k+l,n u k,n 

k, n 

with u -1-s + r s 

k * n k,n k, n k,n+l 

while for some j, r equals the term ratio of the series 
j*n 

we wish to accelerate. 


Constant Splitting Functions 

Since log and arctan are particularly useful, we will honor 
them with yet another example of a new way to 
get an old result; 

n 

- In (1 - x) = x R -x 

n>l n+1 


n + 1/2 2 

arctan x « x R - (_ x ) 

n>0 n + 3/2 

These are just special cases of the hypergeometric function 
with 2nd argument ■ 1, . 

n + a 

F (a, 1; c; z) = R - z , 

n>0 n + o 
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Here m ■> 8 and again, our asymptotic term ratio is z so i f 
we continue to restrict s to be a constant, l/(l-z) is still 
-k,n 

our best guess to minimize u . This follows from setting 

k, n 

s and s to s, r to z, and u to 0 

k, n k,n+l k,n k,n 

in the definition of u 

k, n 

Ue choose j ■ 0 and have 


n + a 


0, n n + c 


1 1 n + a 

1-z 1 - z n + c 



r 

1, n 


n + a 
n + c + 1 


Since this is just r with c+1 replacing c, and c is arbitrary, 
8, n 

we can immediately conclude that 


u 

k, n 

Final Iy. 


k + c - a z 

R - 

k>0 k + c z - l 


n + a 1 

R -, . - 

n>0 n + c 1-z 


c + k - a z 
n + c + k z-1 
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or more symmetrically, 


n z n 2 

- 2 R -z . R - 


n >a n + d z - 1 n>d n + a z - 1 

(with d for c-a), a result usually found by analytic 
methods (See Abramouitz and Stegun, Handbook, of 
Mathematical Functions, Formula IS.3.5). 

(Note that the same transform on this new series yields the 
old one again, so that for every z for which convergence was 
improved, there was the point z/(z-l) for which it was 
degraded, and vice versa.) 

In particular, we have 


x k + 1 x 

arctan x - - R _ _ 

2 k>0 k + 3/2 2 

X + 1 X + 1 

the interval of convergence having expanded from + 1 to 

2 2 2 

± infinity. Setting y . x /(x + 1) w e get 


arcs in y 

2 1/2 
(1 - y ) 


y R - 

k>l k + 1/2 


where 


k 2 •r—\ y 

y . 2^-- 


2 k + 1 


k + 1/2 
k>0 ( ) 

k 


J 

( ) 
k 


( ) 
j-k 


k! (j-k)! 
which, for integer k 

j (j —1) ... (j-k+1) 


k! 

so that - 


2 n 2 k 
( ) ( ) 

n - 1/2 2k k 

( ) , - 

k n k 

( ) 4 













Integrating the arcsin equation: 


2 

2 k + 2 

arcsin y 

2 -Z-f 

y 

k + 1/2 

k>0 

(2 k + 2) { ) 

k 

2 

2 

y 

k 

■ — ■ 

R - 

2 

k>l (k + 1/2) (k + 1) 


2 

y 


2 


y k 2 

- R -y 

2k k>l k + 1/2 

This result, although knoun, is peculiar in the light of 


arcsin y 


■ £ 


k - 1/2 
( ) 
k 


2 k + 1 

y 


2 k + 1 


y 

2 k + 1 


k + 1/2 2 

R -g 

k>0 k + 1 


where i n the coefficients are nearly the reciprocals 
of those in the previous sum. The relationship 
is even more striking if ue write 


y k 2 

arcsin y - — R - y 

2k 1 k + 1/2 
k>- 
2 

which says that if we offset the index of summation by 

2 

arcsin y, we get y arcsin y. 


-1/2 in the formula for 











Other Sp!itting Functions 


Actually, for all but the simplest term ratios, most splitting 

rnmnr 0 2J' const f nt or otherwise, will lead to prohibitively 
comp Iicated u s after just a few k. 
k, n 

So the main problem is to find an s such that u 

k.n k,m 

stays simple. The easiest way 

to do this seems to be to start with r and choose the s 

0, n 

which results in the simplest r you can find. Then repeat 

1, n 

the procedure for r , seeking the s leading to the simplest 
l.n - l, n 

possible r . With some skill and luck, you find a 
2,n - 


sequence s , s , 
0.n l,n 


• • which produces a sequence of u’s 


whose general formula is discernable. 
For example, let us attack 


zeta 2 


1 1 1 

— + — + — 

2 2 2 
1 2 3 


R 

n>l 


(n+1) 


ll pr °'; e convenient to choose j - 1 in the Acceleration 

rormuia, and so define 


l.n 2 

(n + 1) 

certainly can t try s - Lim l/(l-r ), since the asymptotic 
l.n n n 

term ratio is 1, so let us try a very simple function of n 

s «an + b 

1 f n 

with a and b to be determined so as to simplify 
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2 

(1-a) n + (2-a-2b) n + 1-b 

u - - 

1, n 2 

(n + 1) 

A good heuristic is to reduce the degree of u, which we 
can do twice by choosing a = 1, b . 1/2, so that 


l.n 2 

2 (n + 1) 


Then 


2, n 


2 

n 


2 

(n + 2) 


Trying the same form again, 


i - a n + b, 
2, n 


we have 


u 

2, n 


2 

(l-3a) n + (4-4a-4b) n + 4-4b 


2 

(n + 2) 


Analogously, we choose a - 1/3, b - 2/3, so that 


u 

2,n 


4 

2 

3 (n + 2) 


Then 


2 


3, n 2 

(n + 3) 

which suggests 


2 


k, n ' 2 

(n + k) 
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Sticking w i th 

s - a n + b, 
k t n . 


2 2 

(1—a(2k-l)) n + (2k-ak -2bk) n + 


k f n 2 

(n + k) 


which says 

1 3 k - 2 

a - -, b « - 

2 k - 1 2 (2 k - 1) 

3 

k 

u ■ - 

k, n 2 

2 (2 k - 1) (n + k) 

and, sure enough, 

2 


k+1,n 2 

(n + k + 1) 

So finally, 

2 n + 3 k - 2 


k,n 2 (2 k - 1) 

and 


3 

3k. k 

zeta 2 - - R _ 

2 (2 k - 1) k>l 

2 (2 k - 1) 


3 k 1 

. — R -. 

2 k k>l (k + 1/2) 4 


with a convergence rate of 2 bits/term. 


2 

(l-b)k 


2 

(k + 1) 












But this is just a special case of our arcsin formula 
with y - 1/2, so that 


2 

3 8 2 1 pi 

zeta 2 - ~ - arcsin - = -, 

2 2 2 6 . 

(I am indebted to Dr. John U. Urench, Jr. of the Naval Ship 
Research and Development Center for pointing out that this 
zeta 2 formula is derived as a Harkov transformation in 
Section 156 of Knopp’s Theory and Application of Infinite 
Series. Later, however, we will improve this to a 6 bits/term 
formula, which has a better chance of being undiscovered.) 


This example has nicely illustrated the actual practice of the 
method: 

Acceleration Recipe 

1) Write the term ratio.r 

0,n . 

2) Guess a form for a splitting function s , cunningly 

8, n 

including some undetermined coefficients. It probablu should 
be asymptotic to l/(l-r ). 

0, n 

3) Choose these coefficients to siniplifly 


0,n 0,n+l 0,n 

u 

0, n+1 


u 0, n 

0,n 

that it agrees with r and r 

0,n l,n 


u 

k, n+1 

r and determine the 
u k, n 
k, n 

to achieve agreement with r 

k, n 


0, n 


4) Urite r 

1, n 


5) Guess r so 
k, n 


5) Urite r 

k+1, n 


coefficients in s 
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This all boils down to searching for rational solutions to 
the recurrence relations defining u and r , 

k, n k,n . 

and there may be much to gain from exploring the literature of 
nonlinear difference equations, which I have as yet neglected. 


Now let’s try harmonic numbers. 

1 1 1 1 p n 

h - - R 

P 1 2 3 p n>l n+1 

We have not equipped our transform for finite series—it 
would evidently produce a sum of two series o.f orphans, one 
from the left and one from the right. We can, however, use 
the transform to see what the one on the left will be. 

I f we again try 

s - a n + b 
1, n 

wi th 


1, n n + 1 


we find 

n - b + 1 

u » - 

l»n n+1 


and 


2,n n+1 


independent of a and b! This means that we can extract 
an arbitrary series of orphans from the left without 
affecting the term ratio at all.. 

Other choices of splitting functions are equally 
unenlightening. So it wants an infinite series? We 
give it an infinite series. 
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1 1 1 

h » - + - + - + 

P 12 3 

11 1 

- ( — + - + - + ... 

P+1 p+2 p+3 



n>l 



P 

(n + p) 


A slightly less dubious way to get this is to let m approach 
infinity in the cute identity 


m 

P 


n (n + p) 

n«l 

but this seems to rely on p being an integer—a restriction 
which we will profitably waive in order to conquer the di gamma 
(Psi) function while we are at this: 



■ £ 

n«l 


n (n + m) 


h 

P 


£ 


p 


n (n + p) 


Psi(1 + p) + enigma 


uhere enigma is Euler’s constant - .57721SGS49... and is usually 
spelled gamma. We have 

n (n + p) 

(n + 1) (n + p + 1) 

s - a n + b, 

1, n 


1, n 

Using good old 


we get 








PAGE 24 


(1-a) n + (p-2b-a+2) n + (1-bHp+l) 
(n+1) (n + p + 1) 


which suggests 


a « 1 and b 


(p - 1) (p + 1) 

j - - 

1 * n 2 (n + 1) (n + p + 1) 


n (n + p) 


2,n (n + 2) (n + p + 2) 
Let us leap to the conclusion that 
n (n + p) 


k»n (n + k) (n + p + k) 


2 2 

(l-(2k-l)a) n + (p-2kb-(k +jk-p)a+2k) n + (l-b)k(p+k) 


so that if 


(n + 1) (n + p + 1) 


1 3 k + p - 2 

- and b » - 

2 k - 1 4 k - 2 


2 n + 3 k + p - 2 

s - - 

K.n 2 (2 k - 1) 


k (k + p) (k - p) 


^* n 2 (2 k - 1) (n + k) (n + k + p) 


then 













and lo and behold 
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n (n + p) 


k+l,n (n + k + 1) (n + p + k + 1) 

Final Iy 

P n n + p 


P P + 1 n>l n + 1 n + p + 1 


P 

- s R u 

P + 1 k, 1 k>l • k, 1 


P 3 k + p k (k + p) <k - p) 

“—~~ —--—-- R __________ 

p + 1 2 (2 k - 1) k>l 2 (2 k - 1) (k + 1) (k + p + 1) 


P (3 k + p) k - p 1 

- R - 

2 k (k + p) k>l k + 1/2 4 


k - p - 1 

(3 k + p) ( ) 



k>l k (k + p) { ) 

k 


As a consistency check, we can divide through by p, then 
let p approach 0, yielding the zeta 2 equation we got earlier 

Note that if p is an integer > 0, the k-p in the numerator 
of the ratio will terminate this series on the jth term. 

Ue have our finite series back, only it converges faster! 

For instance, 

1111111 7E1 

1 +- + - + - + - + - + - + - m - , 

2 3 4 5 E 7 8 280 


44 49 238 5 4E 13 29 1 

3 15 1E5 12 585 138B 45045 51480 


P 

- R r 

P + 1 n>l 1,n 
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Now 
i t i 
But 
(h 

39 


if we are computing harmonic numbers as exact nationals, 
s easier to add up the first series than the second, 
if we want big ones, this becomes impractical, 
involves 83 digits.) If we resort to decimal approximat 


ons, 


uexan still sum the harmonic series using approximate reciprocals, 
but even this grows impractical as p gets to be a few hundred. 

There is a divergent series involving log p, Euler’s constant, 
and Bernoulli numbers, but even with an abundance of digits of 
these, you cannot extract many more than 2 pi p/ln 10 digits. 

The formula we have just derived seems 

to solve our problem to the tune of 2 bits/term, but on closer 
inspection, we see that the power of p in the numerator 
of the term ratio is one greater than in the denominator, 
so that If p is large, the series will not start to converge 
until k reaches (p - 21/5. J 

To get the most out of this formula, we should notice a degree 
of freedom in the series transform that we have not yet 
exploited namely, that setting n to m in the s and u of 

k,n k,n 

the transformed series corresponds to starting the original 
series at the mth term. (He must also remember to scale both 
series so that the mth term of the original series is of the 
desired size.) Thus we have 


11 1 ■ 

h - - + — + — + ... 

m,p m m+1 m+2 

( 1 1 1 
m+p m+p+1 m+p+2 

- h - h 

m+p-1 m-1 



y;------- 

4 n (n + p) 


n>m 


n>m 
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p (3 k + 2 m + p - 2) k (k + pi (k - p) 

R - 

2 m (m + p) (2 k - 1) k>l 2 (2 k - 1) (k + m) (k + m + p) 

p (3 k + 2 m + p - 2) k (k + p) (k - p) 1 

R --- _ 

2 m (m + p) k>l (k + 1/2) (k + m) (k + m + p) 4 


k + p k - p 
p(3k + 2m + p + l)( )( ) 

k k 


2k k + m k + m + p 
2 m (m + p) (2 k + 1) ( ) ( ) { ) 

k k k 

Now if m IS about as large as p, ue see that this neu term ratio 

IS about 1/8, so that ue have an almost 1 digit/term method 

for summing the reciprocals betueen half a trillion and a trillion, 

for instance. But ue can use this same trick for those betueen a 

quarter trillion and half a trillion, etc., so ue really 

can accurately compute very large harmonic numbers subject to 

a cost factor of log p. 

2 

This example shous that once ue successfully accelerate a series, 
we have a handy way to sum finite but large intervals of that 
series, since the starting value of the summation index is merely 
a parameter in the neu series. In fact, it need not be an integer: 

Ue knou 



k>l z - k pi 


then 

pi z cot piz.-zh -1 
.-z,2z 


■ L 












a 2 bits/term cotangent formula. 
Similarly for Catalan’s constant: 
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1 

2 

7 



h h 

■ 1/A,p 1/2,p 

p->0 IS p 4 p 


2 



k>0 k + 1/4 2k 

( ) ( ) 
k k 


3 

3 (k + 1) 

- - R - 

2 k>0 2 

(k + S/4) (4 k + 2) 


2 

Pi 

8 


Later on, when ue accelerate eta 2, we will get a single 
series formula for Catalan's constant. 


The next example will emphasize a very important point. Ue 
recall that our original goal was to accelerate the particular 
series 

1 1 ' 1 

- : — +-*— +-+ . .. 

1 (1 + p) 2 (2 + p) 3 (3 + p) 
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with 

n (n + p) 

r - -.--- 

l.n (n + 1) (n + p + 1) 

but in the process, ue accelerated the more general form 
n (n + p) 

r - ___ 

k,n (n + k)'(n + p.+ k) 

This means that if ue uish to accelerate some other series 
*ts term ratio is expressible as 


r 

a, n+b 


for some a and b, then we can use the transform ue already 
summing from k - a instead of 1, and setting m (the initial 
of n in the original series) to b+1 instead of 1. Thus if 
uish to accelerate the apparently very different series 


4 


pi 


2 -12 -112 
1 + (■—) + (— _) 

2 2 4 


-113 2 -113 

(-) + (- 

2 4 B 2 4 E 


(Jolley’s Summation of Series, formula #274) 


2 n - 1 2 

ue find the term ratio . (_) 

2 n + 2 

But this is just r with p = 0: 

3/2,n-1/2 


2 

2 n - 1 2 n 

R (-) „ R - 

n>0 2 n + 2 n>-l/2 2 

(n + 3/2) 


3 k - 3 

2 (2 k - 1) 


R 

3 

k>- 

2 


3 

k 


3 

(2 k - 1) 


3 (2 k - 1) 
8 


3 

(2 k + 1) 

R - 

k>l- 2 

32 k ik + 1) 


and 


found, 
va I ue 
ue 


2 

) 
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so 


4 

Pi 


12 


<2 k - 1) k 


z 

k>l 


8 k 


2 k 3 
( ) 
k 


a formula uhich has among its oddities the property of speedina 

2 p ;us^erm ? h Si0Uin9 d0Wn ’ t0 ’ ts eventual convergent ra e°of 
2 b, ts/term. This ,s the only pi series I can recall whose terms 
(and partial sums) are all finite decimal fractions, being integers 
over pouers of 2 (a fact obscured, in this case, by R notation?? 


Many thanks are due Dr. R 
Inc., who first suggested 


chard Pave I Ie of Perception Technology 
transforming the harmonic series. 


Let us now return to 

pi/4 . arctan 1 - 1 - 1/3 + 1/5 - 1/7 + ... 
to see how we can improve on Euler's transform. He have 
2 n + 1 


0, n 


2 n + 3 


l/2 Ch A P veru C s?Lie' SUS0est ! n 9 a fitting function approximating 
1/4. A very simple non constant such function is apparently 

n + a 

s - - 

8. n 2 n + b 

with a and b chosen to simplify 


2 2 2 

4 (b-2a+l) n + (2b +10b-4ab-12a+4) n +3b +5b-4ab-Sa 


0, n 


(2 n + 3) (2 n + b) 12 n t b t 21 

Reducing.. the da ^ e of the numerator as before, we let 

2 

2 ((2 a-2)n+2a -2a- 1) 

U ■ — — — _ 

8,n (2 n + 3) (2 n + 2 a - 1) (2 n + 2 a + 1 ) 

Instead of simply taking the obvious choice of a - 1 
let us proceed directly to ’ 
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0, n+1 


1, n 


0, n 


0, n 


2 n + 1 2 n + 2 a - 1 


(2 a - 2) n + 2 a 


2 n + 5 2 n + 2 a + 3 


(2 a - 2) n + 2 a 


2 a 


to see which choices of a allow at least one of the denominator 
terms to cancel a numerator term. For each pairing of 
a term from the denominator with one from the numerator we 
scale them so that the coefficients of n are equal, then solve 
Mirc^MA a “ h ' ch " akes the whole terms equal. Uhen using 
nALbYnA, the very powerful symbolic math program at HIT 

°" e „ P ?AK,? r T th ' S 0peration succinctly by invoking 
the RESULTANT function, uhose existence and usefulness 
were revealed to me by Professor Hal Abe I son. The resultant of 
the numerator and denominator of the above expression is 

1G 2 2 2 

- 2 a (a - 1) (a - 2) (a - 3) (a + 1) (2 a - 1) (2 a - 3) 

and in general, the resultant of two polynomials 

ih t h ei r coefficients which is 0 when they have a 
common factor or their leading coefficients simultaneously 
vanish. They are painful to compute, but will no doubt come 
back into popularity nou that mach i nes wiI I do them for us. 

The.series in question is so easy to accelerate that any of these 
choices of a WI I I lead to a ratio which is in turn easy to split, 
etc., and in this way we can generate dozens of pi formulas of 
varying convergence rate. Perhaps the most interesting, yet easily 
Then d f ° rmU 3 ' es alon 9 the P a th indicated by choosing a . 3/2. 


l,n 


(2 n + 2) (2 n + 3) 
(2 n + 5) (2 n + 6) 


The simplest guess of r consistent 
k, n 

with k - 0 and 1 is apparently 

(2 n + k + 1) (2 n + k + 2) 
r - __■_ 

k ’ n (2 n + 3 k + 2) (2 n + 3 k + 3) 

which is, in fact, supported by 
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2 n + S k + 3 

s - -- 

k.n 4 n + 8 k + 4 

(which can be found with some pain by determining a and b 
n + a 

in s - - so that r conforms to our conjectured 

K.n 2 n + b k+l.n 

formula for r ). Then 
k, n 


(k + 1) (2 k + 1) (2 n + k + 1) 

u - - 

K.n (2 n + 3 k + 2) (2 n + 3 k + 3) (2 n + 3 k + 4) 

so that 

P' 5 K + 3 (k + 1/2) (k + 1) 2 

■ - R ------ 

4 4 K>0 (k + 4/3) (k + 5/3) 27 


(5 k + 3) (2 k) ! k! 


k+1 

K=^ 2 (3 k + 2) ! 

at better than 1 decimal digit/term. Interestingly, this 
series can be simplified by rerunning the transform just 
one step on this new series, with 

(B k + 8) (k + 1/2) (k+1) 2 

r - -- 

B.K (5 k + 3) (k + 4/3) (k + 5/3) 27 


6 k + 4 

S m - 

B.K 5 k + 3 

happens to be the reciprocal of the old s 

k, 0 


B (k + 1) (k + 2) 
u - ___ 

0 .K 3 (3 k + 4) (S k + 3) 


and 

which 

Then 

so 
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(k + 1/2) (k + 3) 2 

r - -,- 

1-k (k + 5/3) (k + 7/3) 27 

and thus 

P' 1 5 (k + 1/2) (k + 3) 2 

— - 1- R - 

A 24 k>0 (k + 5/3) (k + 7/3) 27 

or 

E S 

--- 

k 3 k + 4 

k>0 2 (2 k + 1) ( ) 

k + 2 

But this formula is just uhat we uould have gotten had we 
accelerated the series 

P' 1111 

— 1 --+-+ - . . . 

4 3 5 7 3 

and then restored the initial 1: which is to say that 

* P u can be simplified further than s R u 

k>1 ka k.B k>0 k,0 

This suggests that when we get an acceleration formula, we 
should check to see if there are any preferable starting 
values of the summation index which differ by an integer 
from the the one initially intended. The brute force way to do 
this is to rewrite 


k,m k>0 k, m 


8,m k>0 


k, m+1 


k ? m 


k, m 


and then take the resultant with respect to k of the numerator 
and denominator of the expression on the right to determine 
the values of m for which it reduces. 

Actually, the fact that there uere two ways to find that 

h,',t P ''!' catlon in the f °rego i ng pi/4 formula was no coincidence, 
but rather a manifestation of a beautiful duality between a 
series and its transform: y oerween a 
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I f 


k, n 


splits the series uhos 


ratio is r to make the series 
k, n 


with ratio r 

n, k 


k+1, n 

- u , then s ■ 1/s 

s k,n n,k k,n 

k,n 


splits 


r back to r . (Note the switching subscripts.) 
n, k k, n 

This can be verified directly from the recurrences defining 
r and u 
k t n k,n 

Thus that extra split by the reciprocal of s that ue performed 

k,0 

on the transformed series was equivalent to the orphaning of 
the first term of the original series via the dual transform. 

Just to make sure the reader understands this duality, an example; 




As we will later show, 
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n + 3 k + 1 

s - --- 

k, n 2n + 4k + l 

spIits 

2 n - 1 


k,n 2 n + 4 k + 1 

to make 

2 k + 1 2 k + 2 

u - - - 

k,n 2 n + 4 k + 1 2 n + 4 k + 3 

so that the transformed series will have ratio 

3 k + n + 4 2 k + 1 2k + 2 

r rn _ __ _ 

n,k 3k+n+l 4k+2n+3 4 k + 2 n + 5 

By duaIity then, 

4 n + 2 k + 1 

s - - 

k,n 3 n + k + 1 


sp I its 


k, n 

to make 


3 n + k + 4 2n + l 2n + 2 

3 n + k + 1 4 n + 2 k + 3 4n + 2k + 5 


u 

k, n 

with the ratio 


n, k 


3 n + k + 2 2 k - 1 

3n + k + l 4 n + 2 k + 3 

2 k - 1 
2 k + 4 n + 1 
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Zb ta 3 


iii i 

zeta 3 - — + — + , V"'' __ 

3 3 3 2-j 3 

123 n>l n 



k>l 


k-1 

(-) 


3 2k 
k ( ) 

k 


the improvement being asymptotically 2 bits/term versus 0 
b i ts/term. To be honest, I first got this formula from a 
Kummer transform (omitted for brevity), so it is probably knoun. 
Ue have 


3 


0, n 3 

(n + 1) 


if we choose s - a n + b, we find that a must be 1/2 to 
0, n 

make the degree of the numerator of u less than the degree of 

0, n 

the denominator. But then we will find that there i s no value of 
b for which r is again 
1, n 

a ratio of cubics, and experience indicates that the larger 

the degree of the term ratio, the harder it is to find a splitting 

function which keeps the degree from increasing yet further. 

The next simplest splitting function asymptotic to n/2 is a 
quadratic divided by a linear polynomial, which will 
contain 3 coefficients to be determined. Now we could begin a 
massive search through alI values of these coefficients for which 
r simplifies, but first, let me suggest two heuristics useful 


for limiting such searches— 
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1) Choose the denominator of the splitting function s to have 
factors in common with the denominator of the term ratio. 

2) Choose the numerator of 1 - s to have factors in common with 
the numerator of the term ratio. 

Uhile these heuristics are not foolproof, a glance at the formula 
defining u will show why they tend to produce simple u’s. 
k, n 

If we use these heuristics to eliminate all three bugger 
factors, we get 

2 

n 

s - 1 +- 

0,n 2 (n + 1) 

which rewards us with 

2 

n (n + 1) 

r - - - 

1* n 2 

(n + 2) (n + 3) 

suggesting 

2 

n (n + k) 

r ■ -'_ 

k.n 2 

(n + k + 1) (n + 2 k + 1) 

But there is no reason to believe that s should be asymptotic 

k,n 

to n/2 for all k. To find out what it should be, we 
temporarily set s to a n and find that a should be 
k, n 

1/(2 k + 2) in order to reduce the degree of u . Now the 

k,n 

factors of the numerator and denominator of r are not all the 
same, so which do we choose for s? We need another heuristic: 

3) Choose the largest factors. 

This dictates 

- 2 

(n + k) 

s = 1 + 

k, n 


2 (k + 1) (n + 2 k + 1) 
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wh i ch 


and s 


l*Jh i ch 


Again 
with i 
term. 


where 
In an 
zeta 3 


Use 


is indeed the winner among the eight candidates. Then 

2 2 

(k + 1) (n + k) 

u - ---'___ 

k, n 2 

(n +. 2 k + 2) (n + k + 1) (n + 2 k + 1) 


i nee s « 5/4, 

k.l 


5 

zeta 3 - - 

2 

4 k 

is the R uay of saying 


- k 

R - 

k>l 2 (2 k + 1) 



k-1 

(-) 


3 2k 
k ( ) 

k 


2 

we notice a similarity to the arcsin formula, this time 
I " 5/2, but with with an extra factor of 1/k in each 
which we can explain by dividing by y and integrating: 


zeta 3 



2 

(asinh y) 

- d y 

y 


asinh is inverse hyperbolic sine. 

effort to rewrite this as the usual integral definition of 
, I tried substituting 

y - sinh t/2 . 


coth t/2 









and 
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fix) dx . / f (x) dx 

Jo 



f(x+a) dx 


and 


zeta(n;x) 


1 


(n - 1) ! 


where 


zeta(nj x) 


k 



k>l k 


n-1 

x t 

-dt 

t 

e - x 


the generating function of the -nth powers. lzeta(n;l) * zeta n 
and zeta(l;x) - - log (1 — x)). Then expanding the numerator 
into 3 terms and using 


2 

1 pi 2 

zeta(2; —) --(In 4>) 

2 15 

4 > 

(see J. Eduards' Treatise on Integral Calculus, p. 2S5+), 
one uinds up with 

2 

5 1 pi 53 

zeta 3 - - zeta (3; —) + — log $-(log®) 

4 2 6 6 

<t> 

where 4 is the golden ratio 


2 

- 1.618... . $ - 1 


asinh 1/2 
e 


Surprisingly, this formula appears in problem XXVII 41 Vii 
of Edwards, credited to the ingenious manipulations of Landen. 
Landen and Ramanujan independently found the analogous relation: 


2 

8 1 2 pi 

zeta 3 - - zeta(3; -) + - log 2 

7 2 84 


1 3 

— (log 2) 
42 
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uhich suggests that there is at least one more of these identities, 

2 

wi th d> instead of 4 or-2. Perhaps this is what was garbled to 
produce the obviously false equation in problem 41 iv. 


Clausen, in Crelle’s Journal, vol. 5, (1830), p38, gives 


P' n 3 n 

zeta 3 - - R - _ _ r _ 

2 1 2 (2 n + 1) 2 n>l 2 (2 n +■ 1) 

8 n n>- 8 n 

2 


Noting that these R’s are integrals of (arcsin y)/y and 


2 

(arcsin y) /y respectively, ue can perform the analogous 
manipulations (starting with y = sin t/2i) to get 


n 



h>0 (3 n + 1) 


3 

5 pi 13 zeta 3 

- + - 

3/2 36 

2 3 


and 


n 



n>0 (3 n + 2) 


3 

5 pi 13 zeta 3 


3/2 36 

2 3 


by separate consideration of the real and imaginary parts. 
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Returning to the zeta 3 acceleration formula, ue note that setting 
n to 1/4 instead of 1 gives us a formula for 


1 1 11 

— + — + — + - + 

3 3 3 3 

1 5 9 13 



so 


1 11 

3 3 3 

1 3 5 


3 

1 zeta 3 pi 

- (zeta 3 -- + -) 

2 8 32 


7 zeta 3 
15 


3 2 

Pi 80 k + 112 k + 41 


54 2 

40 (k + 1) (4 k + 1) 


2 

- 15 (k + 1) 

R - 

k>0 (8 k + 9) (8 k + 13) 
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Splitting Function Addition Formula 

Sometimes it is desirable or even necessary to choose a sequence 
of splitting functions whose formulas alternate between two different 
forms. The most useful special case of this is when 

every other s is identically 1. This has the effect of orphaning the 
first term of each partially accelerated series, replacing n by n + 1 
in -the term ratio. As we shall see, this is useful for increasing 
convergence rates, but more important, it allows us to use 
acceleration formulas which would otherwise involve division by 0. 

Consider again the problem of accelerating zeta 3, with 

3 


0, n 3 

(n + 1) 

Before finding better heuristics, I found, by analogy with.the 
Kummer transform, 

2 
n 

2 (n - 1) 


1 


(n - 1) (n -f 1) 
2 

n (n - 1) 


2 

(n + 1) (n + 2) 

don’t exist, a problem which we might sidestep by 

remembering to sum from n « 2 instead of 1, adding back in the 
initial term when we finish. Blithely proceeding, we wi I I eventually 
discover the pair 

2 

■ n (n - k) 


k, n 2 

(n + 1). (n + k + 1) 


s 

0, n 

g i v i.ng 


u - 

0, n 


and r 

1, n 


But s and u 


k, n 


2 - 


n 


2 (k + 1) (n - k - 1) 
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w i th 


2 

(k + 1) 

u - ---- ( ■ 

k,n (n + k + 1) (n.- k - 1) 

But this says that if we start with n=2, we will get a rero 
denominator in the k=l orphan. So we must translate the 
k«l series by a term as we did for k « 0, In fact, we will 
have to do this for every k. A nice way to think of this 
is that before the application of the 

kth splitting function s , we interject the splitting function 
k, n 

s « 1. Thus, before the k*0 split, r will become 
n 0, n 

3 

(n + 1) 


3 

(n + 2) 


necessitating the replacement of n by n + 1 in s , which 

0, n 

evades the division by 0. By extension,, r will become 

k, n 

2 

n (n + k) 


2 

<n + k + 1) (n + 2 k + 1) 


just before the k + 1st replacement of n by n + 1, 
after which the value of s must be 
k, n 

2 

(n + k + 1) 


2 (k + 1) n 


(i.e. n + k + 1 replaces n). Ue note that the r formula now 
corresponds to the one we got using the heuristics, but the s 
formula does not. This is because the s in the earlier acceleration 
combined the effects of the current s and the s = 1 preceeding it. 
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How could we have combined these two alternating spitting functions 
to get the earlier one, had we not already known it? "Luckily' 1 , 
some straightforward algebra solves the more general problem, of 
finding the splitting function s equivalent to successively 

n 

splitting by p and then q ; 

n n 


Splitting Function Addition Formula 

s - p + q (1 - -p + p r) 
n n n n n+1 n 


In our case, p - 1 and 
n 


q - s 
n k, n 


2 


(n + k + 1) 
2 (k + 1) n 


Using 


our new, 


2 

(n + k) n 


n k,n+k 2 

(n + k + 1) (n + 2 k + 1) 

compound splitting function becomes 


2 

(n + k) 

s ■ 1 -- 

k,n 2 (k + 1) (n + 2 k + 1) 

in accordance with the earlier transform. 

The main point of all this is that once we successfully transform 
a series by finding an s-r pair, we can find any number of new s-r pairs 
by interspersing the unit splitting function using this important 
special case of the addition formula. Ue summarize it under the 
name 







Translation Transform 


s <- 1 + s r 

k * n k,n+k+l k,n+k 


r <- r 
k,n k,n+k 


Rich Schroeppel first proposed the technique of periodically 
translating a series while accelerating it, originally for the 
purpose of improving the convergence rate. Let us use this on 
some accelerations to see how this improves on the improvement. 
First we derive another pi formula from arctan 1: 

• 2 n - 1 


0,n 2 n + 1 

(with n starting at 1), HeuristicaI I y choosing the 
same denominator for the splitting function, while 
requiring it to approach 1/2 for large n, we get the form 


0, n 2 n + 1 

giving 

2 ((2 - 2 a) n + 2 - a) 

U m --- 

0»n (2 n + 1) (2 n + 3) 

suggesting a - 1 so that 

2 n - 1 


1. n 2 n + 5 

suggesting 

2 n - 1 


k,n .2 n + 4 k + 1 

requiring a « 3 k + 1 (in the formula for u not shown) 

k, n 
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so that 

n + 3 k + 1 

s . -.--—■ 

k, n 2 n + 4 k + 1 

giving 

(2 k+ 1) (2 k + 2) 

u - -:_ 

k, n (2 n + 4 k + 1) (2 n + 4 k + 3) 

Then 

P' 3 k + 2 (k + 1/2) (k + 1) 1 

— - —. R .., 

4 3 k>0 (k + 5/4) (k + 7/4) 4 

with a modest 2 bits/term. (The author is mildly embarrassed to later 
discover that this formula is merely the result of pairwise grouping 
the terms of Euler’s transform of the arctan 1 series, and could have 
been computed directly by using p = 1/2 - q in the addition formula.) 

n n 

Using the Translation Transform, we get 

2 n + 2 k - 1 

r ' - - 

k.n 2 n + S k + 1 

n + 4 k + 2 2 n + 2 k - 1 

s = 1 - - - 

k.n 2 n + E k + 1 2 n + 6 k + 3 

These give 

pi 2 (7 k + 6) k + 1/2 k + 1 1 

. 4 15 k>8 k + 7/6 k + 11/6 27 

with a healthy log 27 - 1.43 decimal digits/term. Ue could have 
gotten this pi formula directly by accelerating arcsin 1/2, for which 

(2 n + 1) (2 n + 2 k + 1) 
r - ___ 

k.n 4 (2 n + 2 k + 2) (2 n + 6 k + 3) ’ 


(2 n + 2 k + 1) (2 n + 2 k + 3) 

s ■ 1 ---- 

k.n 3 (2 n + 6 k + 3) (2 n + 6 k + 5) 


and 
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2 

2 (2 k + 1) (k + 1) (2 n + 2 k + 1) 


k - n (n + k + 1) (2 n + G k + 3) (2 n + G k + 5) (2 n + 6 k + 7) 

starting n and k at 0. 

Recalling the fact that offsetting n by 1/2 changes arcsin y 
2 

into (arcsin y)/y, we have 

2 2 
pi 7 k + 5 - (k + 1) 

— . - R - 

18 3 (2 k + 1) k>0 3 (3 k + 4) (3 k + 5) 


2 

9 


E 


k 3 

(-) (7 k + 5) k! 


(2 k + 1) (3 k + 2) ! 


This last formula can also be found by a straightforward acceleration 
and translation transform of eta 2: 













PAGE 48 

(one half times Markov's ceta 2 formula). 

Then, translation transforming; 

n + k n + 2 k 

k,n n + 2k + l n + 3k + l 


s 

k, n 


1 


(n + k) (n + 2 k) (n + 4 k + 3) 

2 (n + 2 k + 1) (n + 3 k + 1) (n + 3 k + 2) 


2 

(k + 1) (n + k) (n + 2 k) 

u - - 

k ' n (n + 2 k + 2) (n + 3 k + 1) In t 3 k t 21 (n + 3 k + 3) 

whereupon s R u gives us 3/2 times the earlier formula 

0,1 k>0 k,l 


2 

for pi /18. Curiously, s factors, so that 
k,5 


2 

pi 1.111 

— 

12 1 4 9 IB 



k 3 

(-) (k + 3) (7 k + 12) (k + 1) ! 


2 

k (2 k + 3) (3 k + 4) ! 


Repeating the translation transform on the pi/4 formula, 
ue get a series with G bits/term, but these terms have an 
unfactorabIe quadratic in their numerators. Noting, however, 
that r is the term ratio for In 2, ue qet 

0,3/2 

In 2 ■ s R u 

k,3/2 k>0 k,3/2 

14 k + 11 (k + 1) (2 k + 1) 

. - R - 

IE k>0 8 (4 k + 5) (4 k + 7) 


3 (14 k + 11) 


k>0 (4 k + 3) (2k 


4 k + 1 k+2 

1) ( ) 4 











Trans I ation 

r 

k, n 

s 

k, n 

u 

k, n 

Here we must 
is initially 

zeta 


providing G 


transforming the zeta 2 series, we get 

2 

(n + k) 


(n + 2 k) 


(2 n + 5 k) (n + k) 


2 (2 k - 1) (n + 2 k) 


3 3 

k (n + k) 


2 2 

2 (2 k - 1) (n + 2 k) (n + 2 k + 1) 


and sum from k - 1 so that r 


that of zeta 2: 


k,n 


21 k - 8 

- R 

8 (2 k - 1) k>l 


21 k 8 


8 (2 k - 1) (2 k + 1) 
3 


8 k>l 3 

8 (2 k + 1) 


1 (21 k -r 13) k! 

8 z—/ 3 

k>0 (2 k + 1)! 


>i ts/term. 











Translation transforming the 4/pi formula, 


2 

(2 n + 2 k - 1) 


k, n 2 

(2 n + 4 k + 2) 


2 

4 n + 10 k + 7 12 n t 2 k - II 

s - 1 +- -- 

k,n S (k + 1) 2 

(2 n + 4 k + 2) 


2 2 
(2 k + 3) (2 n + 2 k - 1) 


k.n 2 2 

8 (k + 1) (2 n + 4 k + 2) (2 n + 4 k + 4) 


Noting that s reduces to 

k,2 

42 k + 47 


32 k + 32 

we get 

3 

4 1 42 k + 5 (2 k + 1) 

— ■ 1 + - +- R - 

Pi 4. 2048 k>l 3 

512 (k + 1) 


5 

4 


2 k - 1 3 

(42 k + 5) ( ) 



k>l 2 


uhich, like its predecessor, is also a series of exact binary 
fractions, its kth term being an integer of about Gk bits 
shifted about 12k bits to the right of the binary point, for 
a net gain of G bits/term. Jack Holloway (of the M. I. T. 
Artificial Intelligence Lab) points out that this enables us to 
compute the millionth digit of 1/pi without computing the 
first half million digits. 
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Translation transforming the pi formula whose ratio approached 
2/27 (where we had set out to beat Euler's transform), 
we have 

(2 n + 3 k + 1) (2 n + 3 k + 2) 


fc.n (2 n + 5 k + 2) (2 n + 5 k + 3) 


(2 n + 3 k + 1) (2 n + 3 k + 2) (2n + 7k + 5) 

s - 1 - —---- 

k.n 2 (2 n + 5 k + 2) (2 n + 5 k + 3) (2 n + 5 k + 4) 


- - (k + 1) 2k + 1 2n + 3k + 1 2n + 3k + 2 2n + 3k + 3 

u - - - - - 

k,n 2n + 5k + 2 2n + Sk + 3 2n + 5k + 4 2n + 5k + 5 2n + 5k + 6 

which promises a term ratio approaching 54/3125. But this 
rapid convergence is marred by the fact that s has an 

k,0 

unf actorab I e cubic numerator. Fortunately, s does not, 

k, 1 

so we can again pull the trick of transforming the series with 
its first term borroued. Thus 

P' (Ilk + 13) (17k + 12) 1 2k + 1 3k + 3 3k + 4 3k + 5 

— - 1 - R - 

4 720 k>0 5 5k + 7 Sk + 8 5k + 9 5k + 11 

or 


Pi 




(Ilk + 13) (17k + 12) 


5k + S k 

k>0 3 (2k +1) (k + 1) (2k + 3) ( ) (-2) 

2k + 3 
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Translation transforming the special cases of the hypergeometric 
function near the beginning of this paper., we get a feu familiar 
identities and a few weirdos. Recalling 


n .+ a n + b 

F (a, b; c; z) » R —_ _ 

n>0 n + c n + 1 

For F(a,l;c;z) we will have 

n + k + a 


k,n n + 2 k + c 

n + k +. a z 

1 -- -- 

n + 2 k + c z - 1 

Then 


s 

k, n 


2 

(z-2)k+(c-a)z-c k +a k + c - a z 

F (a,l;ciz) - - r - - - 

C (z - 1) k>0 2k + c + 1 2k + c + 2 z - 1 


Taking the resultant with respect to k to find simplifying values of 
a and c, we find that the degree of the term ratio reduces in the cases 
c - Z a - 2, 2 a - 1, 2a, 2 a + 1, and 2 a + 2, but only in 
the case c ■ 2 a is the result expressible as an ordinary hypergeometric 


z - z z 

F ta, 1; 2a; z) - -F (a, 1 j a+1/2;-) 

2 Z - 2 4z-4 

In the other cases, the only way to assure the n+1 in the denominator 
of the term ratio requires particular values of zs 
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Fta.lt 2a-2: 


2a~4 

a-3 


(a - 2) 

F(a-2,2; a-1/2; -) 

(a-3)(a-1) 


2 

(a - 1) 


Fta.l; 2a-l; 


2a-3 

a-2 


2 

(a - 3/2) 

F(a-1,2; a+1/2; —-) 

(a-2)(a-1) 


(a - 1) (2a - 1) 


2 

3 (a-1/2) 


Fla,2; a+-; -) 

2a-l 2 a (a-1) 

F(a,1; 2a+l; -) = - 

a (1 - a) (2 a - 1) 

2 

3 a 

F (a,2; a+-; --) 

_ 2a 2 (a-1)(a+1) 

F(a,l; 2>va+2;-} 

a+1 


(1 - a) (a + 1) 
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It is a manipulative pain in the neck to confirm a fact uhich 
is obvious uhen you think about it: 

Translation transforming an s-r pair and its 
dual produces equivalent formulas. 

This is because the orphans split off by the l's are made by 
the s’s in the dual, and vice versa. 


For some reason, discoveries rarely proceed along the most direct 
line. It is only afterwards that one sees how he could have found 
his result much more easily. An example is the next three zeta 3 
formulas, which are derived in seemingly the reverse of logical order. 
This presentation has two advantages over Gaussian sanitization, 
in that besides providing a datum on human problem solving, 
it shous off the power of the splitting function 
addition formula in two unusual applications. 

3 


0, n 3 

(n + 1) 

Employing the heuristics, we try 

3 

In + b) n 

s « 1 +- . 

0, n 3 

a (n + 1) 

a must then be 2 to reduce the degree of u, uhile b = 2 zeros the 
leading coefficient in an ugly quadratic factor of u’s 
numerator. The result is 


3 

(n + 1) (2 n + 5) 

r m - t 

1, n .3 

(n + 3) (2 n + 3) 

Only slightly discouraged by the fact that the ratio is now 
quartic, we should obviously try 

3 

(n + b) (n + 1) 

s - 1 +-:- . 

l.n 3 

a (n + 3) 
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But an exhaustive application of the FACTOR and RESULTANT functions 
of FIACSYNA fails to turn up any promising leads. Since 
zeta 3 formulas are so hard to find, ue desperately try 


3 

(n + b) (n + 1) 

s - 1 +-----, 

l.n 2 

a (n + 3) (2 n + 3) 

Again a must be 2 to reduce the degree of u, and amazingly enough, 
b « 4 gives 


2 

(n + 2) (n + 3) 


2,n 2 ’ 

(n + 4) (n + 5) 


a simplification to something more like the initial ratio. This suggests 
the possibility of a sequence of ratios alternately cubic and quartic, 
generated by an alternating sequence of splitters. If ue succeed, ue 
can use the addition formula to combine pairs of consecutive splitters 
to produce the compound splitter which takes us from one cubic form 
directly to the next. Anticipating this, we relabel the indices so that 
the cubic ratios have consecutive k’s, uhile the quartic ones are 
distinguished by a prime ('). Thus, this last ratio (a cubic) uill be 
renamed r , uhile the previous quartic uill be called r* 

. 1 * n 0,n 

Now we can guess that 


2 

(n + 2 k) (n + 3 k) 


k,n 2 

(n + 3 k + 1) (n + 4 k + 1) 


and we will shortly find that 

2 

(n + 5 k + 2) (n + 2 k) (n + 3 k) 

s - 1 +---- 

k,n 2 

(4 k + 2) (n + 3 k + 1) (n + 4 k + 1) 


gives 
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r ’ 
k, n 


<n + 2 k + 1) (n + 3 k + 1) (2 n + E k + 5) 
2 

(n + 3 k + 3) (n + 4 k + 3) (2 n + S k + 3) 


and then 


s’ 
k, n 


. indeed confirms 


2 

(n + 5 k + 4) (n + 2 k + 1) (n + 3 k + 1) 

2 

(2 k + 2) (2 n + E k + 3) (n + 4 k + 3) 


2 

(n + 2 k + 2) (n + 3 k + 3) 


k+1, n 

(n + 3 k + 4) (n + 4 k + 5) 


Now we combine s and s’ in the addition formula and get a 
staggering expression for the compound splitting function 
involving an unfavorable polynomial of degree six in both k 
and n. Proceeding to compute u anyway, we determine that each 
of these ungainly terms will contribute a generous twelve bitss 

2 

^ 63 k + 70 k + 13 (k + 1) (3 k + 1) (63 k + 50) 

8 (2 k + 1) (3 k + 2) 

132 (3 k + 2) (4 k + 3)"” 

4 

(3 k + 1) (k + 1) 

U - --- 

k.l 2 2 

16 (3 k + 4) (4 k + 3) (4 k + 5) 


4 3 2 

24570k +64161k +62152k +26427k+41S4 
zeta 3 - - 

3456 (2k + 1) (3k + 1) (3k + 2) 


4 3 2 4 2 

(24570k +641Glk +G2152k +26427k+4154) k! (2k+l)l 

2 

48 (2 k + 1). (3 k + 1) (3 k + 2) (4 k + 3) I 


4 

(k + 1) 

R - 

k>0 2 2 

16 (4k+S) (4k+7) 
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This combination of cumbersome expression and rapid convergence 
is reminiscent of too many applications of the translation transform. 
Noting that the translation transform can be inverted, we ask 
what acceleration formula can be translation transformed to get 
this monster. 


Inverse 

Translation Transform 


s - 1 


k,n-k-1 

k, n 

r 


k,n-k-1 

r <- 

r 

k, n 

k, n-k 


These relations follow from substituting n-k for n and n-k-i for n 
in the translation transform, then solving for r and s 

k,n k,n 

In the case in question, s is too large to exhibit, but the result 
k, n 

of the inverse translation transform is a barely tractable 

2 

(n +. k) (n + 2 k) 


k, n 2 

(n + 2 k + 1) (n + 3 k + 1) 


s 

k, n 


then 


2 3 2 2 3 2 

(n+k) (n+2k) (2n +(20k+12)n +(68k +83k+25)n+80k +150k +93k+19) 


4 (2 k + 1) (n + 2 k + 1) (n + 3 k + 1) (n + 3 k + 2) 


(k + 1) (40 k + 29) 

s - 1 +- 

k, 1 2 

3G (3 k + 2) 


4 

(k + 1) (2 k + 1) 


k,l 2 2 

' 9 (2 k + 3) (3 k + 2) (3 k + 4) 


and 
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so that 

.2 4 

364 k + 501 k + 173 (k + 1) 

zeta 3 - - R --- 

144 (2 k + 1) k>0 2 2 

9 (3 k + 4) (3 k + 5) 

2 6 
__ (3E4 k + 501 k + 173) k! 


k>0 36 (2 k + 1) (3 k + 2) ! 


which is somewhat nicer and stiil provides nearly 3 digits/term. 
Reapplying the inverse transform, 


2 . 

n (n + k) 

r - -- 

k, n 2 

(n + k + 1) (n + 2 k + 1) 

(k + 1) (n + 3k + 2) n + 3k + l 

s * — — --—*——— +--- 

k.n 2 2 (2 k + 1) 

4 (n + 2 k + 1) 


so that 


9 


30 k + 19 (k+1) 

zeta 3 - - R - 


IB (k+1) (2 k + 1) k>0 2 

4 (2 k + 3) 


-£ 


38 k - 11 


3 2 k 2 

4 (2 k - 1) k ( ) 

k 


which, at 4 bits/term, is similar to the very first zeta 3 formula 
we got, but twice as fast. It is, however, quite different from 
and computationally superior to grouping terms of the slower series 
pairwise, which produce an unfactorable cubic numerator. 
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In general, grouping terms pairwise tends to double the 
degree of the term ratio while doubling the digits/term, 
so if you intend to evaluate an R expressed series sequentially, 
the doubled computational effort per term will countervale 
the greater convergence rate. Also, it greatly simplifies machine 
computat i on of these series when the integers comprising the 
numerators and denominators in the R expression do not exceed 
the computer's register size. 

On the other hand, Rich Schroeppel has noted that n 
iterations of pairwise term grouping will, 

n 

for large enough n, evaluate 2 terms of a sum much more rapidly 
than sequential ev-aluation, assuming the Fast Fourier (n log n) 
multiplication algorithm. 

Combining terms pairwise is neatly done in R notation; 

... (b + r (b + r ( ,.. 

2k ' 2k 2k+l 2k+l 

(b + r b + r r (... 

2k 2k 2k+l 2k 2k+l 

or 


- (b + r b ) R r r 

2k 2k 2k+l k>0 2k 2k+l 


b R r 
k k>0 k 
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Pitfalls 

Inverse transI ation. transforming arcsin 1/2, 
2 n + 1 2 n - 2 k + 1 


k,n 8 n + 8 2 n + 4 k + 3 


8 n 2 n + 1 


3 2 n - 2 k.- 1 2 n +' 4 k + 3 


2 k + 1 4 k + 4 4 k + 4 


2 n - 2 k - 1 2 n + 4 k + 3 2 n + 4 k + 5 


2 k + 1 4 k + 4 4 k + 4 

r - - -- -'_ 

n,k 2 k - 2 n + 3 4 k + 2 n + 5 4k+2n'+7 ’ 


Noting that this simplifies for n » 1, and that arcsin 1/2 begins 
11 

- (1 + — ( 1 + ... 

2 24 

so that a is 1/48, 

0,1 

pi 1 1 4 k + 4 4 k + 4 

62 30 k>0 A k + 7 4 k + 9 

But this is absurd because it says that pi > 3 + 1/5, since the 
R expression is > 1. Uhat went- wrong? Recall that in the 
derivation of the transform, the kth partial sum of the orphans 
comprising the new series differs from the original sum by 
the kth partially accelerated series 


s 

k, n 

g i ves 

u 

k, n 



n>l 


a R r = a u u 

k, 1 n>l k,n 0,1 0,1 1,1 


u R r 

k,1 n>l k,n 


which we assumed went to zero with increasing k. 
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Now as k approaches infinity, r approaches 

k, n 


2 n + 1 
16 (n + 1) 

which is the term ratio for 


- 1/2 

(1 + 1 / 8 ) 


by the binomial theorem, so 


2 n + 1 1 

R - . 1- 

n>0 16 (n + 1) 1G 


3/2 

2 n + 1 2 

R - . - 

n>l 1G (n + 1) 3 


Thus the R expression in our error term approaches 
3/2 

2 

1G ( 1-) 

3 


while the u's which are its coefficient form a convergent 
infinite product! He can, in fact, evaluate this product as 
the discrepancy between the two sides of the false equation. 
To do this, ue must find out uhat 


1 4n + 4 4n + 4 

30 n>0 4n+7 4 n t 3 


really equals. Writing 

4 n 4 n + 4 k 


k,n 4 n + 4 k + 3 4n + 4k + S 
n + k 

s « 1 +- 

k,n k + 1 

Then 

n + k 4 k t 3 4 k t 5 


k.n k + 1 4 n t 4 k t 3 4 n + 4 k + 5 


(whose product over k, we note with relief, "diverges to zero" 
n > -1) and 


for 
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4 k + 3 4 k t 5 

1. k 4 k + 7 4 k + 9 

But this is just the sum 

111111 

- (- --■ + - +- ... 

2 3 5 7 9 11 

uith its terms grouped pairwise. Using 

1/2 

2 x 

arctan - 

2 3 5 7 9 

1 - X X X X X 


1/2 3 5 7 9 

2 


evaluated at 1, ue get 


1 

R 

4 

n + 4 

4 

n + 4 

30 

n>0 

4 

n + 7 

4 

CD 

+ 

C 


Final Iy, then. 


1 Pi 

- (- 1 ) 

2 3/2 

2 


2 k t 1 4 k + 4 4 k + 4 

n —-— - - 

k>0 2 k - 1 4 k + 5 4 k + 7 


1/2 

3 2 pi 


8 


This relation can be verified by multiplying by 


sin pi/4 
p i /4 


n 

k>0 


4 k + 3 
4 k + 4 


4 k + 5 
4 k t 4 


to leave 


2k + l 4 k t 3 
k>0 2 k - 1 4 k + 7 
whose kth partial product telescopes to 
2 k + 1 3 

- 1 4 k + 7 

which approaches -3/2. More generally, we have the 
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Corrected Series Transform 


R r , “ s R . u + nu R r 

n>m j.n k.m k>j k.m k>j k.m n>m n 

where r is the limit of r as k grows infinite, 
n k, n 

In our particular case, II u converges for all m, giving 

k, m 

2 k + 1 4k + 4 4 k + 4 

n - - - 

k>j 2 k - 2 i + 1 4 k + 2 m + 3 4 k + 2 m + S 

as a sum of two sums divided tag a third sum. I wish I had time 
to further investigate this connection between sums and products. 

Hoping for another product formula hy dualizing the transform 
that started all this, ue are very unpleasantly surprised. 


2n + l 4n + 4 4 n + 4 

r - - —_ _ 

k.n 2 n - 2 k + 3 4n + 2k + 5 4n + 2k + 7 


3 2 n - 2 k + 1 4 n + 2 k + 3 

s - - - - - 

k.n 8 k 2 k + 1 


2 k + 3 2 k - 2 n - 1 

u - - - # 

k.n 8k 2 k + 4 n + 5 

The product of the u "diverges to 0", leading to a contradiction 
k, n 

of the fact that the two series differ by a finite product. 

The problem here is very serious: each individual split is invalid! 
The quadratic growth of s is too fast for the input series, 
k, n 

which converges like a reciprocal quadratic for all k. 

It was most unnerving and enlightening of Eugene Sal amin 
of the Draper Laboratories to point out that any series can be 
converted to any other series in a single split, since s can 

n+1 
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always be chosen to satisfy the recurrence 

1 - s + r s = anything . 
n n n+1 


Since each of our series transforms involves applying an infinite 
number of splits, how can we be sure ue aren’t perpetrating the 
wiIdest absurdities? Our salvation lies in zippers. To see 
what actually happens to our series, imagine a zipper 
traveling along, splitting it in two. Right behind it, another 
zipper is zipping it back together, but offset so that one tooth 
(the orphan) is left hanging off the front and there is a 
corresponding extra tooth on one side between the zippers. 

If the series were finite, this extra tooth would become the 
orphan on the right when the zipper reached the end. In symbols, 
when the zippers reach the pth term, 



B-l 

a s + \ ' (1 - s ) a + s a 

0 0 n n n+1 n+1 

n=0 


+ (1 - s ) a + \ > a . 

P P n 

n>p 

If the original series converges, then the last 
sum will vanish as p grows, but this isn’t necesariIy true of 
that extra tooth, (1 - s ) a . Thus, besides worrying about 
P P 


the product over k of the u , we must make sure that the 

k, n' 

limit as p grows large of (1 - s ) a approaches something 

k, p k, p 

(preferably 0) which we can some over k. At last, then, the 
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Summing the Reciprocal Fibonacci Numbers 
By an involved process, I found 
k 

(-) f 


k, n f 

n+2k+l 


k 

f + (-) f 

n+Ak+2 n+2k+l 


k,n (2 f + f ) f 

2k 2k+l n+2k+l 

n+k+1 

(-) f f 

2k+l 2k+2 


k, n (2 f + f ) f f 

2k 2k+l n+2k+l n+2k+2 

where f , f , f , are the fibonacci numbers 0, 1, 1, ... 

0 1 2 


so that 


1 

Sr ■ 

n>l n 


k 

f + (-) f k 

4k+3 2k+2 (-) 

- R - 

f f k>0 2 f + f 

2k+l 2k+2 2k+2 2k+3 


k(k-l) 


k 2 

(f + (-) f ) (-) 
4k+3 2k+2 



k>0 2k+l 2k+2 1 3 2k-l 


n In 

where g ■ 2 f + f - 4 +(--), the nth Lucas number, 

n n-1 n 4 

This is an extremely rapidly convergent series, yielding about 

2 

k /5 digits for k terms. 
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Similar manipulations yield 


Z 


1 


sinh nx 


z 


sinh (2 k + 2) x -sinh (4 k + 3 
sinh(2k+l)x sinh(2k+2)x (2-2cosh x) (2-2cosh 3x) 


k>0 

flore genera I I y , if 


n n 

f - x + c y 
n 

so that 

f - (x + y) f - x y f 

n+1 n n-1 


then ue have 


k k 
x y f 


k,n f 

n+2k+l 


k+1 


f 

n+4k+2 

(x y) 

f 

n+2k+l 

k k+1 

k+1 

k 

(x - y ) 

(x 

y ) f 

n+2k+l 

n+k 

2k+l 

2k+l 2k+2 

c (x y) 

(x 

y ) tx 

k k+1 

k+1 

k 

(x - y 

) (x 

- y ) f . 

n+2k+l. 


X 


... (2-2cosh(2k+l)x) 


2k+2 

- y ) 


n+2k+2 
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so that 


f 

m+4k+2 


(x - 1) 


Er 

n>m n 


k+1 

(xy) f 


m+2k+l 

m+k 

C (xy) { 

2k+l 2k+l 

x - y 

f k>3 

k+l k+2 

k+2 k+l 

i m+1 

(x - y 

) lx - y 


2 k+2 2 k+2 

(x - y ) 


f f 

m+2k+2 m+2k+3 


again with the. convergence rate of a theta function. 

The author is perplexed by the fact that terms in the right side 

j j+i j+1 j 

of this equation bI on up when x = y or x = y for some 

nonnegative integer j, while the left side does not. In fact, it 
is hard to see how equality could hold when x is fixed at 2, for 
example, and y is nearly 1. Numerical experiments, however, 
indicate that the equation indeed holds, the large term being 
rapidly eroded by later ones. This indicates that when such an 
offending ] exists and is positive, the above expression, with 
R instead of R , must be 0. 
k> j k>0 


n 

If we really do wish to sum l/(x - 1), (a so called Lambert series 
when x « 1/z), we must resort to sneaky tricks: 


2 

E — 

n>l x - 1 


(Ixl > 1) 


-Z 


1 1 


--f--. 

1/4 n n 1/4 n n 

(x ) - (-1) (-x ) - (-1) 


This holds because, on the right, term 2 m + 1 cancels 
term 4 m + 2, leaving only the terms where n » 4 m. 

These latter two sums avoid the division by zero when fed to the 
transform. 
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In other cases ue can avoid the problem by using the identity 


n 

c y k 

( - ) 

n 

x 


n n 

x - c y 


n 

c y k+1 
{-) 

k n 

k y n x 

c (-) . - 

k+1 n. n 

x x - c y 


repeatedly to establish the Rummer transform 


1 1 

n n • 1 k 1 y k 

n>0 x - c y k>0 (-) - - (-) 

c x c x 


(lx! > ]y|) . 
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One Last Zeta 3 Formula 

Knopp, Chapter VIII exercise 128a, points out a relation 
between zeta 3 and the series with ratio 

3 

n 


3 

(n + 3) 

Tasting blood, we propose 

3 
n 

r - --- 

k,n 3 

In t 2 k t 1) 

as a candidate ratio in yet another zeta 3 transform. 
After guessing the form 

4 3 2 

an + bn + cn + dn + e 

s - 1 -- 

k,n 3 

(n + 2 k + 1) 

we find (with just a little help from nACSYnA) 

1 2 k + 1 

a - -, b « - 

6 k + 2 3 k + 1 


3 

(2 k + 1) 

c - - 

2 (3 k + 1) (3 k + 2) 

3 

(k + 1) (2 k + 1) 

d - 

2 (3 k + 1) (3 k + 2) 

2 3 

(k + 1) (2 k + 1) 

e - -.- _ 

3 (3 k + 1) (3 k + 2) 
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Despite these composite coefficients, the author has so far 
fai led to find any nicely factorable expression or sum of 
two expressions for this splitting function. Putting this 
function in some intelligible form could shed much light on how and 
when to look for such things. As it stands, the gross disparity 
in elegance between the ratio and its splitter seems to suggest 
that we should seek a new formulation of the transform process. 

In any case., taking partial fractions with respect to n. 


3 2 2 

n + 4k + 1 (2k + 1) (3n + (15k + 9) n + 20k + 25k + 8) 

s » --- 

k,n 2 (3k + 1) 3 

E (n + 2 k + 1) (3 k + 1) (3 k + 2) 

Then 

2 

5G k + 80 k + 29 

zeta 3 - - R 

3 k>0 

24 (2 k + 1) (k + 1) 

In the previous exercise, Knopp indicates that 


1 2 

V' ■---- - 10 - Pi 

Z_/ 3 3 

n>l n (n + 1) 

This is 1/8 (- the first term) of what we would get by substituting 
k + 1/2 for k in the last R expression for zeta 3, since 


3 

n 

r - -. 

1/2,n 3 

(n + 2) 


2 

(k + 1) 


3 (3 k + 4) (3 k + 5) 


Euler's Constant 

In the American Mathematical Monthly, vol. 76 #3, 

MarG9 p273, I. Gerst and H. F. Sandham give 

*1122223 3 4 

gamma - _ + _ 

2345 678 15 16 

the numerator increasing by one each time the denominator reaches 
a power of 2. He can write this as the double sum 
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1111111 
2 3 4 5 6 7 8 

11111 
4 5 E 7 8 

1 1 1 

6 3 10 


which is the series for - log 2 started at the 2 th term, 
summed for n.> 0. Recall the application of the translation 
transform to 


2 n - i 

r - - 

k,n 2 n + 4 k + 1 

(the "embarrassing" formula) to get pi/4, since for k - 0 this is 
the term ratio for arctan 1. Replacing n by n + 1/2, and redefining 

n 

r - -- , 

k,n n + 2 k + 1 


we have, for k - 0, the ratio we want. Replacing n by n + 1/2 
in the corresponding splitting function and u, 

2 n + G k + 3 

s * - 

k,n 4 (n + 2 k + 1) 


(2 k + 1) (2 k + 2) 

u ■ --- 

k,n 4 (n + 2 k + 1) (n + 2 k + 2) 


which says 
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n 



n>m 


2m + 6k + 3 (2 k + 1) (2 k + 2) 

- R ----- 

4 m (m +1) k>0 4 (2 k + m + 2) (2 k + m + 3) 


- (m-1) ! 

(for m even). Thus, 


E 


(2 m + E k + 3) (2 k)! 
k+1 


k>0 4 (m + 2 k + 1)! 

to get Eulers constant, we must 


sum the R expression with m - 2 , for j > 0. 


j+1 


gamma 


EE 

j>0 k>0 


6 k + 3 


2 + 2 k + 1 

) 

2 k 


( 


j+2k+2 j 
2 (2 


1 ) 


This prospect grous 
cheerful when we notice that after a feu j, the R expressions will 
converge with great rapidity, returning 2 j + 2 bits/term. This 
is because in the term ratio, there are tuo m's in the denominator 
and none in the numerator. Thus b bits of gamma will require about 
(b-j)/(2j+2) terms of the jth series. This says that the total 
number of terms of all series will be about 


b b-1 b-2 

- +-+- 

2 4.6 


0 

2b 


- (In b - 1) 
2 
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Zeta 2 Revisited 


n (n + K) 

r = --- 

k,n 2 

(n + k + 1) 


n + 2 k + 1 


k,n k + 1 


g i ves 

(k + 1) (n + k) 

u - -- 

k, n 2 

(n + k + 1) 

Now s ■ 2 and 
k, 1 


2 

(k + 1) 

u = -- 

k, 1 2 

(k + 2) 

which is a strange way to prove zeta 2=2 eta 2 . 

Trans fat ion transforming, 

10 k - 3 k 

zeta 2 - - R -_ 

4 k (2 k - 1) k>l 2 (2 k + 1) 

This series has similar convergence to our first (Harkov) one 

3 . k 

zeta 2 ■= - R_ - 

2 k k>l .2 (2 k + 1) 

except that the former one alternates. 
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The zeta 2 term ratio 

2 

n 


2 

(n + 1) 


is a speciai case of 
n + p n + q 
n + r n + s 

Choosing the splitting function a n + b will produce a new ratio of 
n + p n + q Q(a,b;n+1) 


n + r+1 n + s + 1 Q(a,b;n) 

Where Q is a polynomial quadratic in n and linear in a and b. 
Now a and b can be chosen to make Q divisible by any two 
of the quantities n + p , n + q , n + r , or n + s , or 
they can be chosen to zero one or both leading coefficients 
of Q. Thus for our new ratio we may have any one of the 
quantitites 

n + p n + q n + p n + q + 1 

n + r + 1 n + s + 1 n + r + 1 n + s + 1 

n + p n + q 

n + r + 1 n + s-1 

n + q n + p + 1 n + q + 1 

- f or - -, 

n+s n + r + l n + s + 1 

or any of the above with p and q swapped or r and s swapped. 
Also, from the addition formula, we learn that p «= 1 and 

n 

q - - 1/r combine to produce 8, the identity split, 
n n 

and thus the splitting function 
n + r - 1 n + s-1 


n + p-1 n + q-1 


n + p 
n + r 


n + p n + q 
n + r + 1 n + s 


will perform a reverse translation (the opposite of s 
producing 


1) 
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n + p-1 n + q-1 
n + r - 1 n + s-1 

By combining the above operations, we see that it is possible 
to independently change any of the four terms of our original 
term ratio by any integer. Thus, any series whose ratio is of 
this form, with p, q, r, and s all integers, must evaluate to 

2 

c pi + d 

ui ^. c , and d rational, since the splitting process only adds and 
multiplies by rationaI quant i t i es. Among these rational quantities 
is 0/0, which figures into those splits uhich would result 
in cancellation to a linear term ratio. For example, to get 

n + p n + q 


n + r + 1 n + s-1 
we need the splitting function 
s-l-r n + s-1 


s-l-p s-l-q 

Thus, if both denominator terms exceed both numerator terms, 
there, is no way to pass through any ratio where a 
denominator term is as large as a numerator term. In this 


2 

way we.are prevented from proving that pi is rational, 
since it is easy to show that any series with ratio 


n + p 


n + r 


telescopes to.a rational sum when with r is rational and 
r - p is an integer. As a corollary, we observe that we can 
achieve such a cancellation if a denominator term initially exceeds 
one in the numerator, thus any series uith ratio 


n + p n + q 
n + r n + s 

sums to.a rational if p is rational, r - p is an integer, and 
q - s is a nonnegative integer. 


A part i cu I ar |.y neat special case of the above 
is the transform: 
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r 


k, n 


2 

n 


(n + c) (n + k) 


n +. c - 1 

s - -- 

k» n k + c - 1 


2 

k 

u « - 

k* n (k + c - 1) (n + k) 

which says 

2 2 
n ni + c - 1 k 

R -- - - R - 

n>m (n + c) (n + j) j + c - 1 k>j (k + c) (k + m) 

i. e. 


i s 


(j + C - X) R 
n>m 


2 

n 


(n + c) (n + j) 


a symmetric function in c, j, and m. 


With ail of these relations, there must be a neat irrationality 

2 

of pi proof in here somewhere. 
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For those readers who have found the treatment so far a bit on the 
Heaviside, warning: we are approaching a singularity—we will now 
attempt to accelerate the "series" 

e E (1) 0! - 1! + 2! - 3! + ... 

1 


gotten from iteratively integrating by parts e dt/t from 1 
to infinity. 

Ue have 

r - - n . 

1, n 

This is a special case of 

n (n + c) (n + d + 1 ) 

k 


k»n (n + a + k) (n + d ) 

k 

with d » c « a and z » -1. 

1 

- 1 

s - - 

k.n z (n + d ) 
k 

splits this to yield 

(1 + z (d + k + a - c)) n + (k + a) (1 + z d ) 


k » n (n + d ) (n + a + k) z 

k 

which formally satisfies the recurrences if 
(k + a) (1 + z d ) 

k 

d » - 

k+1 1 + z (d + k + a - c) 

k 


This says that 
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- 1 

R r - R v 

n > m j.n z (m + cl ) k>j k, m 

J 

where 

1 + z (d + k + a - c) 
k 


k. n (n + a + k) z 


l+z(2k+2a-c-l) (k + a - 1) (k + a - c - 1) 

(n + a + k) z '(n + a + k - 1) (n + a + k) v 

k-1 

(using the recurrence relation for d to eliminate it). 

k 

Expressed in standard notation, 


E a a y^b 

j» n j*m Z—4 k, m 

k> j 


with r 


k, n 


a b 

k.n+l k+1, n 

-- and v - -- 

a k t n b 

k.n k,n 


so that b - 

k+1, m 

1+z 12 k + 2 a - c - II 
(m + a + k) z 


k + a-1- k + a- c-1 

b - --- -b 

k.m m + a + k m + a + k - 1 k-l,m 


w i th b 


■ s and b ■ s v . 

j * m i * m j +-1, m j, m j, m 


In the special case of 0! - 1! + 2! - 

ue have j-m-1, c-a-d - arbitrary, and z - -1. 

1 
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For reasons unknoun to the author, this last series seems to 
oscillate with ever decreasing frequency and amplitude about 
the "correct" answer .536347..., gotten from the integral. 

The convergence is slow: summing through b gives 

76 

.533233... and through b gives .5335.... . 

100 

Larger values of c tend to reduce the frequency of early 
oscillations and the amplitude of later ones, slightly improving 
the long term convergence. 
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Numerical Experiments 

The reason that ue choose s approximately « 1/(1 - r ) 

k.n k.n 

and not exact Iy is that r usually becomes an explosively 
k, n 

complicated expression after a feu k. But there is nothing 
preventing us from numerically grinding out a few orphans 
to see how fast the series would go, could ue find its 
formula. Ue try good old 

pi - 4 (1 - 1/3 + 1/5 - 1/7 + . .. 

and get the following data. (The orphans are consecutively 
numbered, uith the cumulative sum on the next line. The correct 
portion of each sum is underlined. Sums and orphans have a I I been 


S3 

scaled up by a factor of 10 and truncated to integers.) 

0 300000000000000000033B00033C00G00000BG000000380303C3BBB00G8000G00B0000 
30030300003000B0BB00B030BGG0B0B0030003000000000000300G00G0GG3G000000O0 

1 1388888883888888388S883888888S88888888S888888888888883828838888888883 
31388888888S88888S888888888888S88SS888S8S88888883838883S8SS88888888888 

2 26721.6184977729594 54382236447748437629958499321706260091949G7146498 
31415S1050738ESG184834327112533GG3732G51884738821059514838883856035378 

3 313522276894413966379925585183206S033402S578222228694221244 785286 
3141592402966355626230965105092182053302218765399281743592305100820664 

4 249396611247391121684367925111189166634236262217053865571729697 
3141592652362966873622086789460107164431385399635543960651170672550361 

5 1231734176348781444565983771714828110936413086526108111308129 
3141592653594701049970868234026090936206213510571957047177278783858490 

6 - 7800184639717679772350514693761088428831202445735635954831 
3141592653586900865331150554253740421512452422143075844731543147903659 

7 2891758706224834123600203816344000323984340648544949096724 
31415926535897926240373753883773436253287964224670601S5330088097000383 

8 320785467783247360695814560164132363740064579366477048 
3141592653589793544822843171624701321143356586599423925444667463477431 

3 - - 307978293062687771788265106518483534771642225750468795 

3141592653589793236844550108936929532878250068115889153802441713008636 
J 1767408805137966711559705978848431072884617127544505 

3141592653589793238611958914074896244437356046964380226687058840553141 

11 - 241276143030542884996759163025924061191175058234964 
3141592653589793238370682771044353359441196883938456165495883782318177 

12 91705799864633482047373735623914069796610114571307 
3141592653589793238462388570908986841483570619562370235292493896889484 

13 21817387880486815563912123128461308427248586867 
314159265358979323846241038829S86732S3041395316854936966G8321145476351 

14 - 25843147094975436508273013732023811464818370 
3141592653589793238462410362453720233328703023412484964577109680657981 

15 23298B224886S99277781293S778645368S3277459031521 
3141592653589793238462643342678607132606484314230349501445387139689502 


10 
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10 301176019786152588257236371345S7.S39645SS1667 

3141592553589793238462643372796209111221743140813356633825026785241169 
^■ 7 £43452815262727517168061437934656142244005 

314159265358979323846264337363S66192648447E65718201813396368°927485174 
18 - 47883177665521571424556495866729058605936 
314159^.653589793238462643373591778748818943085757461637096953876879268 

19 9688659760934822448413 03821465383865878 3632 
3141592653589793238462643383280438509753771534170439851755792535662870 

20 - 221213928865692447844914858845300630 
3141592653589793238462643383280438288539842668478052005340933693362190 

21 _, ,, - 1808993349637383331178873274388527 

3141592653589793238462643383280438286730849318840671675662860415973663 

22 337363528495513808425916164736 
31415926535897932384E26433832SB4382S67311866823691671894 704S6332138459 

23 - 335556307196828079374140488069782438074 
31415926535897932384626433832795027304229893542891 

24 260780474970340433891581727267086 
31415926535897932384626433832795027336847703292595 

20 ,, rr-rr-r 7763798411391533256336209705386 

31415926535897932384626433832795027306925341276709 
20 - 2523054640525121915223658804 

31415926535897932384626433832795027306925316076162 

27 226172366098368594381506560540833915 
31415926535897932384626433832795029568648377059848 

28 - 1491958619247912513790035935615 
31415926535897932384626433832795029568634057473656 

29 ,, rr.-r.r-~r- 30633096936074657867266587 

31415926535897932384626433832795029568634057779987 
30 „ - 72673167398180031715041434344190950 

314159265358979323846264338327950288419 
31 „ - 1638666569936777925212081330108 

31415926535897332384626433832795028841 
32 , 9816226177585909871383193496275 

314159265358979323846264338327950288419 

33 214159255358979323845264338327950288419 23099889185933100428308 

^314159265358979323845264338327950288410 ^02222673517598341 

35 7305315247817 

314153255358979323846264338327950288419 

3B 3J41E92B53589793238.E2643383279502884:3n6| ,6S ‘ 557326 " S563JU32K1 ” 58 

37 3141592653ES979323S46Z6433832795028S41371633 " 24400703417335193 
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We note several things: 

1) The first.six terms are positive, after which the signs 
are apparently random. 

2) The series converges faster and faster until the 
first negative term, which provides only one 
additional correct digit. 

3) After this, the terms grow erratically smaller until 
term 15, which has nearly 1000 times the magnitude of term 
14 and yet provides four new correct digits. 

4) Terms 18, 22, 26, 29, 31, 33, 34, 35, and 37 have 
the urong sign, i.e. they actually make the sum worse. 

Term 31 actually undoes a good digit. 

5) Term 36 is seventeen orders of magnitude ]_arger than 

term 35 and provides five new digits, while the sizable term 27 
overcompensates and contributes less than one. 

6) The general trend after term 6 seems to be ever slower 
convergence. 

Little wonder that we couldn’t find a formula—I dare 
you to prove it converges at all! 

I am indebted to Richard P. Howell of the HIT RLE PDP-1 
for confirming that these results are not merely the PDP-10 
bignurn routines gone haywire. 

Now in the early 1950’s, Daniel Shanks did a Ph. D. thesis 
(Journal of Mathematics and Physics, 1955, first article) 
involving primarily numeric methods of sequence extrapolation. 

Among them was a transform which he denoted e , 

1 

which, in effect, replaces each but the first term with 
the limit of the exponential sequence defined by that term’s 
predecessor, itself, and its successor — the first member of the 
kth sequence produced by this process becoming the kth member 
of the transformed sequence. On page 5 Shanks gives a 
numerical example for which he too chose the 4 arc tan 1 
series (actually the sequence of partial 

sums). He computed only six transform terms to eight decimal 
places, since he was technologically constrained to electro.- 
mechanical computation. Rich Schroeppel made the remarkable 
observation that if Shanks had started his sequence with 
a zeroth partial sum of 0, his transformed sequence would have, 
been identical to the screwy sequence of partial sums 
that we got from splitting by l/(l-r), and indeed, extending 
the computation of his example (without the initial 0) 
exposes the same erratic convergence. Although it is doubtful 
that he anticipa ted this, Shanks provides some valuable insight 
into this phenomenon with his section on spectra of sequences. 
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n i see Iianeous Accelerations 

In exercise 129c of chapter VIII, Knopp gives a relation- 
completely equivalent to one step of a transform for eta 2 and 
Catalan’s constant, but he seems to miss the significance of it, 
presumably because he uas more interested in identifying series 
than in evaluating them. Rephrased in terms of splitting functions, 

2 


k, n 2 

(n + 2 k) 

i k 2n + 6k + l 

k,n 2 2 2 

(n + 2 k) 

3 

k (2 k + 1) 

u - --- 

k.n 2 2 

2 (n + 2 k) (n + 2 k + 1) 


Inverse translation transforming the "embarrassing" pi/4 
formulas 

r 

k, n 

1 (4 k + 1) (4 k + 5) (4 k + 1) (2 k + 1) 4 k + 1 

s - - +---+- 

k - n 2 1G (n + 3 k + 1) 8 (n + 3 k + 2) 1G <n - k - 1) 


n - k 

n + 3 k + 1 


1 k + 1 2k + l 4k + l 4k + 3 
u - - - - - - - 

k - n 2 n-k-1 n + 3 k + 1 n + 3k + 2 n + 3k + 3 


so i f n starts at 1/2 and k starts at 0, 

Pi 11 k + 8 8 k + 1 

— . - R - -- 

4 15 k>0 3 2 k + 1 


4 k + 1 
G. k + 7 


4 k + 3 
G k + 11 


■ z 

k>0 


k 

(XI k + 8) 16 k! (3 k) ! (4 k) ! 


3 (8 k + 5) (2 k) ! (6 k + 1) ! 


which, with a term ratio approaching 16/27,■has the least 
finite convergence rate of any (rational) pi series I know. 
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Even worse is the 0 bits/term one gets by purposely choosing 
sp I i 11ing functions which don* t make u tend to zero with, n: 

n + k n-k + 2 


k,n n + k + 2 n-k + 1 

2 n + k 

1- - 

4 k + 3 n-k + 1 

4 k - 1 n + k n - k 

u - -- -• _ 

k.n 4k + 3 n + k + 2 n-k + 1* 

Choosing n - -1/2, k = 0 (i.e. starting the arctan 1 series one 

term too early), 

^ Pi 4 k + 5 2 k + 1 4 k- - 1 

4 4 k + 3 k>0 2 k + 3 4 k + 3 


k 

(4 k + 5) (-1) 

Xj (2 k + 1) (4 k - 1) (4 k + 3) 
k>0 

which is no better then if we had taken terms of the original 
series pairwise. Grouping pairwise, however, would have 
destroyed the a Iternating sign. Instead, this series 
arises directly from arctan 1, 

2 k' + 1 


k 2 k + 3 

with the single application of the splitting function 

2 k + 1 


k . 4 k - 1 
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Retransforming Euler’s transform of arctan X, 


0,n 2 n + 1 


I found 


2 

n + k 2 n + 4 k + G k + 3 

r - -' _ 

k,n 2 n + 4 k + '3 2 

2 n + 4 k + G k + 1 

4 n + G k + 3 

s - - 

k,n 2n + 4k + 3 


u 

k, n 


2 k - 1 n + k 

2 n + 4 k + 3 2 n + 4 k + 5 


2 

2 n + 4 k + 14 k + ’ll 


2 

2 n + 4 k + G k + 1 


which is- very unlike the other accelerations in this paper 
due to the quadratic occurrence of k in r and u. 

Translation transforming the foregoing, then initializing 
n to 1 and k to 8, 

2 

pi 8 (14 k + 33 k + 20) 2 k + 2 2 k - 1 


2 105 k>0 3 6k+ll Ek+13 

Noting the similarity between this formula and our ear Iier 

-Pi 2 (7 k + G) k + 1/2 k + 1 1 

4 IS k>8 k + 7/6 k + 11/6 27 

we suspect a single application of some splitting function relates 
them. To find it we divide the first formula by 2 and rewrite the 
second formula so that they both have the same expression to the 
right of the R. 


2 (7 k + G) 

2 

k + 1 2 k + I 

k + 2 

G k 

+ 13 

2 k - 1 

15 k>0 

3 

Sk+7 Gk+11 

k + 2 

G k 

+ 13 

2 k - 1 


2 (7 k + 6) 2 k - 1 6 k + 7 

is. IT" "~T~* 


1 2 k + 2 2 k - 1 

+ 1 k>0 3 6 k + 11 G k + 13 
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b b 

k+1 k 

(Thrice using the rule R -r . — R r .) 

k>0 b k b k>0 k 
k 0 

Now ue have tuo series of the form 

b R r and c R r . 
k k>0 k k k>0 k 

Performing a modified split on the first of these, 

(b-t+t) R r - t R r + (b - t ) R r 

k k k k>0 k k k>0 k k k k>0 k 

-t+rt R r + (b - t ) R r 

0 0 k+1 k>0 k+1 k k k>0 k 


-t+rt R r + (b - t ) R r 

0 k k+1 k>0 k k k k>0 k 


- t +(b - t +r t ) R r 

0 k k k k+1 k>0 k 

which ue must make equal to c R r with the right choice of t . 

k k>0 k k 

This means t =0 and 

t - b + c ' 
k k k 

t . -- . 

k+1 r 

k 

In this case, 

2 

4 (14 k + 33 k + 20) 

b - - 

k 105 

2 (7 k + G) (2 k - 1) (S k + 7) 

c - - 

- k 105 (k + 1) 


2 k + 2 2 k - 1 

k 3 Gk +11 Sk+ 13 ' 
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Factoring the first few numerical values of t , ue can easily guess 

k 

2 k 6k + 5 G k + 7 

t - - - - 

k k + 1 5 7 

uhich is confimed by direct substitution into the recurrence. 

Noting 


(b-t+r t )R r 
k k k k+1 k>0 k 


t t r b b 

k k+1 k k+1 k+1 

b (1 - — +-} R - p 

0 b b b k>0 b k 

k k+1- k k 

we see that t /b is the ordinary splitting function 
k k 

for the series with ratio b r / b . 

k+1 k k 

Thus 


3 k (G k + 5) (G k + 7) 


k 2 

2 (k+1) (14 k + 33 k + 20) 


splits 


P'l 

4 

to get 


pi - 


2 

4 (14 k + 33 k + 20) 

-- p 

105' k>0 


2 k + 2 2 k - 1 

3 G k .+ 11 G k + 13 


2 (7 k + G) 

-- R 

15 k>0 


2 k + 1 2 k + 1 

3 Gk+7 Gk+11 


eliminating that ugly quadratic. Perhaps these bothersome 
polynomials which arise so often in our faster accelerations are 
similarly dispensable. Perhaps there is even a way to do this 
to series whose convergence rate has been doubled several times 
by pairwise term grouping. 
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Negat i ng both sides of the recurrence for t , 

k 


(- t ) 
k+1 


(- t ) - c + b 
k k k 


which says that if t takes b to c , then -t takes c back to b . 

k k k k k k 

Thus 

3 k E k + 5 
k 7k + £2k-l 

carries the second series back to the first. 


Another example.of alternating splitting functions: 
First a preliminary split of 


n 2 n + 1 

by 

4 n + 4 

s ■ - 

n 2 n + 3 

to get 

n 2 n - 3 
0. n 2 n + 7 2 n - 5 

Then 

n + k 2 n - 4 k - 3 


k.n 2 n - 4 k - 5 2 n + 8 k + 7 

i s spIit by 

4 n + 4 k + 4 

p - - 

k.n 2 n + 8 k + 7 

to get 








n + k 


2 n - 4 k - 1 
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k,n 2 n - 4 k - 3 2n+8k+ll 
whichis in turn spilt by 

4 n + IB k + 11 ' 

q - - 

k,n 2 n + 8 k + 11 

to get r . p and q combine in the addition formula to 

k+1 

2 n - 4 k - 3 Gk + 5 6k + 9 

s - p + q - - --- 

k »n k,n 2 n - 4 k - 5 2 n + 8 k + 7 2n + 8k + 9 

which gives a u containing five linear terms over five linear terms 
k, n 

which approaches -27/512. s is quartic, and for at least the small 

k, n 

positive integers, won’t factor. The resulting pi formulas are dull. 
Mildly interesting is the fact that translation transforming 
yields a formula both faster and slightly simpler, but not enough 
to make it a good pi formula. (UnfactorabIe cubic on the left of 
R, four linear factors over four linear factors on the right, 
ratio approaches -27/3125.) 

Somewhat different pi formulas result if the preliminary splitting 
function s is merely 2, which gives 
n 

n 

r - - - r’ 

. 0,n 2 n + 3 -l,n+l 

so that q and p can be combined with the addition 

k-l,n+l k, n+1 

formula to get a different splitting function. 

In general, the arguments to the addition formula my be 
interchanged by making the first argument a preliminary split 
on the first iteration. Unless one of them is identically 1 
(the translation transform), the result!ng formu!as wiI I 
differ significantly, but converge at the same rate. 







